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ABSTRACT 

The  conversion  of  electromagnetic  to  acoustic  energy  is 
a  broad  and  complicated  subject.   Some  topics  can  be  found 
in  the  literature  but  its  thorough  treatment  is  rather 
complex. 

Engineering  approaches  to  the  problem  are  however  possi- 
ble through  useful  simplified  models.   Simplified  treatments 
of  the  three  major  mechanisms  leading  to  the  production  of 
mechanical  forces  on  materials  supporting  an  electromagnetic 
field  are  to  be  presented.   These  mechanisms  are  referred 
to  as  radiation  pressure,  transient  surface  heating  and 
stimulated  Brillouin  scattering. 

A  discussion  of  the  experimental  apparatus  and  procedures 
used  in  an  attempt  to  experimentally  verify  some  of  the 
results  given  in  the  theoretical  formulation,  is  also  made. 
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I.   INTRODUCTION 

A.   NATURE  OF  THE  PROBLEM 

The  generation  of  mechanical  forces  on  materials 
supporting  an  electromagnetic  field  was  first  predicted 
several  decades  ago  although  its  detailed  theoretical 
treatment  had  been  done  much  more  recently.   The  literature 
concerning  this  phenomenon  is  not  plentiful  and  is  rather 
complex.   Since  the  efficiency  of  the  energy  conversion  is 
extremely  low,  the  experimental  verification  of  the  theor- 
etical predictions  was  only  possible  after  high  power 
electromagnetic  sources  became  available.   Radiation  from 
powerful  radars  as  well  as  from  high  power  pulsed  lasers  are 
the  most  suitable  sources  of  electromagnetic  energy  for  the 
study  of  the  several  mechanisms  involved  in  this  conversion 
process. 

The  possible  interpretations  of  the  mechanical  forces 
produced  in  a  material  medium  subjected  to  interaction  with 
an  electromagnetic  field,  may  be  divided  into  three  major 
types.   That  is,  the  forces  may  appear  due  to  one  or  more 
of  the  following  mechanisms: 

(i)   radiation  pressure, 
(ii)   transient  surface  heating, 
(iii)   stimulated  Brillouin  scattering. 

It  is  a  familiar  fact  that  electromagnetic  waves  carry 
energy  through  the  medium  where  they  propagate  in  accordance 
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with  Poynting's  vector.   Less  familiar,  however,  is  the 
fact  that  electromagnetic  waves  may  also  transport  linear 
momentum . 

Originally,  the  term  radiation  pressure  was  used  to 
designate  the  pressure  exerted  on  an  object  by  shining 
light  on  it.   This  was  predicted  theoretically  by  Maxwell 
about  1870  but  the  first  measurement  of  radiation  pressure 
was  made  latter  in  1901-1903.   According  to  Maxwell's  pre- 
dictions, if  a  parallel  beam  of  light  carrying  the  energy 
U  falls  on  an  object  for  a  time  t  and  is  entirely  absorbed, 
the  magnitude  of  the  momentum  "p*  delivered  to  it  is  given 
by  p  =  U/c,  where  c  is  the  speed  of  light.   The  direction 
of  p  is  the  direction  of  the  incident  beam.   When  the  energy 
U  is  entirely  reflected,  the  magnitude  of  the  momentum  de- 
livered is  twice  that  given  above,  or  p  =  2U/c.   On  the  other 
hand,  if  the  energy  U  is  partly  reflected  and  partly  absorbed, 
the  delivered  momentum  will  lie  between  U/c  and  2U/c. 

The  resultant  forces  are  so  small  in  relation  to  forces 
of  the  daily  experience  that  they  cannot  be  ordinarily 
noticed. 

Nowadays  radiation  pressure  has  a  much  broader  signifi- 
cance.  It  is  interpreted  as  the  mechanism  responsible  for 
the  production  of  forces  in  a  medium  due  to  a  rate  of  change 
of  electromagnetic  momentum  associated  with  moving  energy 
as  its  equivalent  mass.   This  rate  of  change  of  electro- 
magnetic momentum  is  the  consequence  of  energy  gradients 


established  in  the  medium  due  to  conversion  of  energy  to 
heat  or  to  progressive  changes  in  the  storage  energy;  it 
can  also  be  the  consequence  of  the  existance  of  velocity 
gradients  as  it  is  the  case  when  part  of  the  incident 
energy  is  reflected  back  at  the  interface  between  two 
different  media. 

The  radiation  pressure  forces  are  in  general  much  smaller 
than  the  forces  produced  through  the  transient  heating  of 
the  surface  of  a  medium. 

The  generation  of  acoustic  waves  utilizing  the  thermal 
effects  accompanying  the  absorption  of  electromagnetic 
radiation  by  an  elastic  medium  was  first  experimentally 
verified  in  1881  and  since  then  just  a  few  theoretical  and 
experimental  works  have  been  developed.   The  subject  seems 
to  be  receiving  now  more  attention  probably  due  to  the  fact 
that  the  recent  advances  in  high  power  laser  technology 
gives  to  this  process  the  possibility  of  concrete  and  useful 
practical  applications. 

When  a  pulse  of  electromagnetic  radiation  passes  through 
an  absorbing  medium,  its  intensity  is  exponentially  reduced 
producing  a  rise  in  the  temperature  of  the  medium.   The 
thermal  shock  occuring  as  a  consequence  of  the  intense  sur- 
face heating  produces  a  thermal  expansion  originating  an 
acoustic  wave  which  propagates  away  from  the  heated  surface. 
That  is,  a  simple  thermoacoustic  source  is  established  in 
the  medium. 
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It  is  possible  to  extend  this  concept  of  a  simple 
thermoacoustic  source  to  the  one  of  a  thermoacoustic 
array  which  may  be  produced  in  the  medium  when  it  is  sub- 
jected to  a  train  of  periodic  pulses  of  high  energy.   In 
this  case,  a  highly  directional  acoustic  wave  is  launched 
from  the  heated  region  and  propagates  in  a  direction  normal 
to  the  electromagnetic  beam. 

As  suspected,  the  analysis  of  this  thermodynamic  con- 
version process  involves  the  knowledge  of  theoretical  prin- 
ciples related  to  fields  like  solid  and  fluid  mechanics 
and  thermodynamics . 

So  far,  the  mechanisms  leading  to  the  conversion  of 
electromagnetic  to  acoustic  energy  were  described  by  an 
interaction  process  between  an  electromagnetic  wave  and  the 
matter  supporting  it.   Under  certain  conditions,  sound  may 
also  result  when  two  intersecting  laser  beams  of  very  high 
energy  interact  with  a  medium.   If  the  angle  of  interaction 
of  the  two  beams  is  the  appropriate,  optically  induced 
generation  and  amplification  of  sound  will  be  verified.   This 
is  a  particular  case  of  the  scattering  phenomenon  and  is 
called  stimulated  Brillouin  scattering.   Note  that  this  is 
not  the  same  mechanism  as  the  transient  surface  heating. 
While  the  latter  utilizes  the  thermalization  of  modulated 
electromagnetic  energy,  the  former  is  based  on  the  mechanics 
of  travelling  wave  parametric  interactions  where  the  resul- 
tant acoustic  wave  has  a  frequency  equal  to  the  difference 
in  the  frequencies  of  the  two  interacting  light  beams . 
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B.   BACKGROUND  CONCEPTS 

A  more  detailed  theoretical  analysis  of  the  major 
mechanisms  responsible  for  the  conversion  process  will  be 
presented  next.   But  before  that,  it  seems  to  be  worthwhile 
to  make  a  brief  but  hopefully  clarifying  comment  about  some 
unfamiliar  terms  belonging  to  fields  like  tensor  calculus, 
mechanics  and  thermodynamics.   The  terms  to  be  introduced 
here  are  mentioned  in  the  theoretical  formulation  to  be 
presented. 

Tensor  analysis  —  this  is  a  rather  advanced  subject  in 
mathematics  and  so  just  a  few  basic  concepts  will  be 
mentioned. 

In  a  three  dimensional  system  consider  two  vectors  A 

-*-  ■** 

and  B  such  that  each  rectangular  component  of  vector  B  is 

a  linear  function  of  the  components  of  vector  A,  that  is 


Bl  -  TUA1  +  T12A2  +  A13A3  ' 


B2  =  T21AX  +  T22A2  +  T23A3  , 
B3  =  T   Aj  *    T32A2  +  T33A3  . 


If  the  association  of  the  components  of  vector  B  with 
those  of  A  is  to  be  maintained  for  any  rotation  of  the 
system  coordinates,  then  the  coefficients  T. .  have  to  trans- 
form  in  a  specific  manner.   By  definition,  a  tensor  of  rank 
two  is  a  linear  transformation  of  the  components  of  a  vector 
A  into  the  components  of  a  vector  B  which  is  invariant  to 
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rotations  of  the  coordinate  system.   The  nine  components 

T.  .  of  the  linear  transformation  are  invariant  scalars 
ID 

and  are  called  the  tensor  components. 

The  divergence  of  a  tensor  of  second  rank  is  a  vector, 
or  tensor  of  first  rank,  and  the  divergence  of  a  vector 
is  an  invariant  scalar,  or  tensor  of  zero  rank.   The  concept 
of  a  tensor  can,  of  course,  be  generalized  to  include  ten- 
sors of  arbitrary  rank. 

Stress  -  suppose  that  a  cut  perpendicular  to  the  x  axis 

is  made  in  a  body  in  equilibrium.   In  order  to  maintain  the 

equilibrium,  an  internal  force  would  have  to  act  in  the 

exposed  face.   By  definition,  the  stress  at  that  point  is 

the  intensity  of  such  a  force,  that  is,  the  force  per  unit 

area.   If  AF  is  the  resultant  contribution  of  these  internal 

forces  on  an  element  of  area  AA,  then  the  components  of  AF 

along  the  three  coordinate  axis,  AF  ,  AF   and  AF  ,  are 
3  x    y        z 

called  stress  components  and  are  defined  by 


AF 

ax  =  lim  -JF 
X    AA+0  AA 


AF  AF 

T     =  lim   — r-£  ,     T     =  lim   -rr- 

xy    AA_0  AA       xz    AA_Q  AA 


a   is  the  normal  stress  and  t   ,  x    are  the  shear  stresses 
x  xy   xz 

and  from  its  definitions  it  is  clear  that  normal  and  shear 
stresses  are,  respectively,  the  intensity  of  force  perpen- 
dicular and  parallel  to  the  plane  of  the  area  element.   Normal 
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stresses  are  termed  tensile  if  they  pull  away  from  the  cut; 
those  pushing  against  the  face  are  called  compressive 
stresses. 

Obviously,  any  deformation  of  a  body  due  to  external 
forces  will  give  rise  to  internal  forces,  or  stresses.   The 
stresses  may  then  be  thought  of  as  entities  describing  the 
state  of  the  internal  action  of  a  body. 

Imagine  now,  as  an  extension  of  this  concept,  that  a 

cube  of  infinitesimal  dimensions  is  isolated  from  a  body 

and  has  its  faces  perpendicular  to  the  coordinate  axis. 

The  stesses  exposed  by  the  cube  are  then  the  three  normal 

stresses  a    ,  a      and  a      and  the  six  shear  stresses  t   , 
x   y      z  xy 

t   /  t  ,    x  ,    t    and  t   .   It  can  be  shown  [Ref.  1]  that 
yx   yz    zy    zx      xz 

t    =  t   ,  t    =  t    and  T    =  T   . 
xy    yx   xz    zx      yz    zy 

The  stress  components  may  be  arranged  into  a  matrix 
format  representing  the  state  of  stress  at  a  point: 


yx 


zx 


xy 


zy 


xz 


yz 


This  is  called  a  second  rank  stress  tensor. 

Strain  —  the  concept  of  strain  is  related  to  the 

deformations  of  a  body.   The  normal  strains  in  the  axis 

directions,  e  ,  e  and  e  ,  are  defined  as  the  change  in 
x   y      z 

length  per  unit  length  of  a  line  segment  in  the  direction 


14 


under  consideration.   It  follows  from  this  definition  that 
strain  is  a  dimensionless  quantity.   Shear  strains  Y- • 
may  be  also  defined;  they  represent  changes  in  angles 
between  line  segments.   The  subscripts  in  y . .    have  similar 
meanings  to  the  subscripts  for  shear  stresses. 

The  strain  components  can  also  be  assembled  into  a 
second  rank  strain  tensor. 

It  should  be  noted  that  stress  and  strain  are  closely 
related  quantities. 

Elasticity  —  is  the  property  of  a  material  which  enables 
it  to  return  to  its  original  size  and  shape  after  the  cause 
of  a  deformation  has  been  removed. 

Bulk  modulus  of  elasticity  'B'  —  also  called  modulus  of 
dilatation,  is  defined  as  the  ratio  of  the  change  in  pressure 
on  a  body  Ap,  to  the  resulting  fractional  change  in  volume 
-AV/V,  or  B  =  -VAp/AV. 

Specific  heat  'C  —  this  is  the  heat  capacity  per  unit 
mass  of  a  body;  heat  capacity  is  the  ratio  of  the  heat  AQ 
supplied  to  a  body  to  its  corresponding  temperature  rise 
AT.   So,  C  =  AQ/mAT. 

Thermal  conductivity  'K'  —  is  a  constant  of  proportion- 
ality relating  the  time  rate  of  heat  transfer  dQ/dt  across 
a  certain  area  A  and  the  corresponding  temperature  gradient 
dT/dx,  or  dQ/dt  =  -KA  dT/dx. 

Mechanical  equivalent  of  heat  'J1  —  since  heat  is  just 
another  form  of  energy,  any  energy  unit  could  be  a  heat 
unit.   Hence,  the  mechanical  energy  equivalent  of  heat 
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energy,  that  is,  the  number  of  joules  equivalent  to  one 
calorie  is  called  mechanical  equivalent  of  heat, 
J  =  4.185  J/cal. 

Coefficient  of  linear  thermal  expansion  'a1  —  the 
change  in  any  linear  dimension  of  a  body  is  called  a  linear 
expansion.   This  coefficient  is  then  a  constant  of  propor- 
tionality relating  the  temperature  change  AT  of  a  body, 
its  original  linear  dimension  I   and  the  corresponding  change 
in  this  dimension  Lt,    or  A£  =  a£AT. 

Adiabatic  process  —  is  a  process  taking  place  in  such 
a  way  that  no  heat  flows  into  or  out  of  the  system.   In 
practice,  since  the  flow  of  heat  is  somewhat  slow,  any 
process  can  be  made  approximately  adiabatic  if  it  is  per- 
formed quickly  enough. 
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II.   THEORY 

A.   RADIATION  PRESSURE 

1.  Elastic  Stress  Tensor 

The  study  and  description  of  a  certain  category  of 
physical  phenomena  like  for  example,  volume  deformations, 
cannot  be  done  correctly  by  the  use  of  the  scalar  or  vector 
concepts.   The  adequate  treatment  of  this  type  of  problems 
requires  the  use  of  the  tensor  concept;  tensor  calculus, 
a  somewhat  unfamiliar  tool,  plays  then  an  important  role 
in  the  analysis  of  such  problems. 

When  dealing  with  mechanical  forces  on  a  material 
medium,  it  is  necessary  to  introduce  an  entity  called  elas- 
tic stress  tensor.   This  is  a  tensor  of  second  rank  and  so 
it  has  nine  components,  the  stresses,  representing  the 
forces  exerted  on  unit  elements  of  area.   It  turns  out  that 
the  divergence  of  this  tensor  is  a  vector,  or  tensor  of 
first  rank,  giving  the  force  per  unit  volume  caused  by  the 
elastic  stresses. 

2.  Electromagnetic  Stress  Tensor 

An  electromagnetic  field  is  now  considered;  as 
usual  in  electromagnetic  theory,  its  definition  is  based 
upon  the  four  vectors  E,  B,  D  and  H  satisfying  Maxwell's 
equations.   E  and  H  are,  respectively,  the  electric  and 
magnetic  field  vectors.   The  electric  and  magnetic  flux 
density  vectors,  D  and  B,  are  related  to  E  and  H  by  the 
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constitutive  equations  D  =  eE  and  B  =  yH  where  £  is  the 
permittivity  and  y  is  the  permeability  of  the  medium  under 
consideration . 

The  existence  of  an  electromagnetic  field  in  a 
material  medium  leads  to  the  definition  of  the  so  called 
electromagnetic  stress  tensor;  its  nine  components  T   ,  T   , 
T   ,  T   ,  etc.   The  electromagnetic  stresses  are  defined 
[Ref.  2]  by  the  following  relations 


T     =   |  (E2  -  E2  -  E2)  +  %    (H2  -  H2  -  H2) 
xx     2    x    y     z     2    x    y    z 


similarly  for  T   ,  T 

■*  yy        zz 


T  =T  =£EE+yHH 

xy  yx  x  y  x  y 


/ 


similarly  for  T   ,  T  (1) 

2  yz    zx 


The  divergence  of  this  tensor,  div  T,  is  a  vector 
quantity  whose  components  in  cartesian  coordinates  are 
obtained  according  to  the  rule 


8T      3T      3T 
(div  T)    =  -^  +   -JEZ.  +   ^i£2 
x      3x     3y     dz 


similarly  for  (div  T)   and  (div  T) 

(2) 
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Using  Maxwell's  equations  it  can  be  shown  [Ref.  3] 
that 


div  T   =   pE  +  JxB  +  e|~  (E"xB)  (3) 


where  p  and  J  are,  respectively,  volume  charge  density  and 
current  density. 

As  already  mentioned,  the  divergence  of  the  elastic 
stress  tensor  is  interpreted  as  a  force  per  unit  volume;  this 
important  result  can  be  extended  to  the  case  of  the  electro- 
magnetic stress  tensor  so  that  its  divergence  may  be  inter- 
preted as  an  alectromagnetic  force  per  unit  volume. 
3.   Electromagnetic  Momentum 

It  is  apparent  that  the  terms  pE  and  JxB  in  the 
right-hand  side  of  eq.  (3)  have  dimensions  of  a  force  per 
unit  volume.   For  a  medium  with  p  and  J  throughout  a 
volume  V,  these  quantities  represent  then  the  forces  exerted 
on  unit  volume  elements  of  charge  and  current,  that  is,  they 
are  the  well  known  electric  and  magnetic  forces  exerted  on 
charges  and  currents .   This  means  that  a  portion  of  the  total 
force  acting  on  the  unit  volume  appears  as  a  force  on  the 
matter  within  that  volume.   Therefore,  the  remaining  force 
exerted  on  the  unit  volume  should  be  associated  with  the 
term  e  8/3t  (Exb)  . 

As  the  dimensions  of  B  =  uH  are  MQ   T   and  those 
of  5  =  £E*  are  QL   ,  the  quantity  e  (Exb)  is  dimensionally 
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-2   -1  -1      -2  -1       -1  -3 

QL  MQ  T   =  ML  T   =  MLT  L   or,  a  momentum  per  unit 

volume  since  the  dimensions  of  momentum,  p  =  mv,  are  MLT 
As  a  consequence,  the  term  e(ExB)  is  interpreted  as  the 
electromagnetic  momentum  density.   Now  it  is  clear  that  the 
remaining  force  exerted  on  the  unit  volume  by  the  external 
field  acts  by  increasing  the  electromagnetic  momentum  within 
that  volume. 

The  precise  manner  in  which  electromagnetic  momentum 
density  should  be  defined  has  found  some  differences  of 
opinion.   The  existence  of  a  duality  of  treatment  of  such 
concept  results  mainly  from  the  two  different  ways  in  which 
matter  can  be  visualized.   On  the  one  hand,  a  macroscopic 
presentation  of  a  medium  may  be  adopted  when  dealing  with 
electromagnetic  propagation  in  that  medium;  on  the  other, 
a  microscopic  picture  describing  matter  as  an  arrangement  of 
dipoles,  doublets  and  charges  surrounded  by  vacuum,  exists 
associated  with  the  concept  of  energy  whose  value  depends  on 
the  relative  velocity,  but  which  always  moves,  regardless 
of  the  medium,  with  the  velocity  of  light  in  vacuum. 

According  to  the  latter  interpretation,  the  electro- 
magnetic momentum  density  is  defined  [Ref.  4 ]  by  y  e  (ExH) , 


where  c  =  l//y  e    is  the  velocity  of  light  in  free  space. 

To  summarize,  it  can  be  stated  that  associated  with 
an  electromagnetic  field,  there  is  a  momentum  distributed 
with  a  density  p  which  is  given,  according  to  the  specific 
interpretation,  by  one  of  the  two  possible  expressions 
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p   =   yQ£o(gxH)  (4) 


p   =   ye(ExH)  (5) 

Throughout  this  work  the  latter  interpretation  will  be 
mostly  used. 

4 .   Forces  in  the  Electromagnetic  Field 

The  divergence  theorem,   a  useful  mathematical 
relation  in  electromagnetic  theory, 

£      D-dll   =    /   V-D  dv  (6) 

S  V 

states  that  the  integral  of  the  normal  component  of  any 
vector  field  over  a  closed  surface  S,  is  equal  to  the 
integral  of  the  divergence  of  this  vector  field  throughout 
the  volume  enclosed  by  the  closed  surface. 

This  vector  divergence  theorem  has  a  tensor  analogue 
stating  that  the  integral  of  the  divergence  of  a  tensor 
throughout  a  volume  V  is  equal  to  the  integral  of  its 
normal  component  over  the  surface  S  enclosed  by  V. 

Integrating  eq.  (3)  over  the  volume  V  and  applying 
this  theorem,  the  result  would  be 


&        T.dS   =   /   pE  dV  +   /   (JxB)dV  +  £  %-r      f    (E*B)  dV 

*   V 
S  V  V  (7) 
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Consider  for  the  moment  a  stationary  distribution 

of  charge  and  current;  the  fields  are  independent  of  time 

and  so  the  last  term  on  the  right  of  eq.  (7)  is  zero.   The 
other  two  terms 

/  pE  dV  =  f     ,      /   CJxB)  dV  =  F  (8) 

e  m 

V  V 

represent  the  net  mechanical  forces  acting  on  the  volume 
distributions  of  charge  and  current,  where  p  and  J  are  the 
charge  and  current  densities  within  the  closed  surface  S. 

The  total  force  transmitted  by  the  electromagnetic 
field  across  the  surface  S  or,  the  resultant  force  exerted 
on  the  charged  matter  within  S  is 


Fe  +  $m  =   /   (pli  +  ^Tx£)  dV  (9) 


It  can  be  shown  [Ref.  3]  that 


F^  =  /   PE"  dV  =  /    [eC£-n)2  -  |  E2n]  dS       (10) 


e 


V  S 


t       -  /   (jxS)  dV  =   /   [  i  (S-n)S  -  ~  B2n]  dS   (11) 
m  y  2y 

V  s 


where  n  is  the  outward  unit  normal  at  a  point  in  the  surface 

enclosing  V,  E2  =  E2  +  E2  +  E2  =  |e|2   and 

B2  =  B2  +  B2  +  B2  =  |B|2  . 
x    y    z    '  ' 
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By  Newton's  second  law 

?e  +  K     m     ft  5m  <12> 

P*m  is  the  total  linear  momentum  of  the  charged  matter  within  S 

Consider  now  that  the  volume  V  contains  no  charges 
or  currents.   If  the  field  is  stationary,  the  net  force 
across  S  is  zero.   But  if  the  field  is  variable,  the  third 
term  in  the  right-hand  side  of  eq.  (7)  is  not  zero  and  so 
this  equation  becomes 


6   T-dS   =   ~    /   £(ExI)  dV  (13) 

S  V 


The  electromagnetic  momentum  concept  shows  up 
now  to  help  in  understanding  this  apparent  action  of  a  force 
arising  in  neutral  matter  immersed  in  an  electromagnetic 
field.   As  stated  before, e (ExB)  =  p  is  the  density  of  the 
electromagnetic  momentum  associated  with  the  electromagnetic 
field.   The  integral 

/   p  dV  =  Pe  (14) 

V 

is  the  total  momentum  of  the  field  contained  within  V  and 
then,  as  eq.  (13)  says,  the  force  transmitted  across  S  goes 
to  increase  the  electromagnetic  momentum  of  the  field 
within  V. 


23 


Consider  at  last  a  system  composed  of  charges  and 
variable  fields  within  a  bounded  region.   The  conservation 
of  momentum  theorem  is,  therefore,  according  to  eq.  (7) 
expressed  by 


4r   (Pm  +  ?J   =  t     T-dS  (15) 

dt    m    e 


It  follows  from  eqs .  (9) ,  (10)  and  (11)  that 


f     =   /   [e(E-n)E  +  y(H-n)H  -  j    (eE2  +  yH2)n]  dS     (16) 


This  is  the  resultant  force  on  the  charge,  current  and  matter 
within  the  surface  S . 

Equation  (16)  is  a  valid  relation  for  both  stationary 
and  dynamic  situations  and  the  validity  of  its  application 
to  the  case  of  variable  fields  can  be  justified  if  basic 
principles  of  the  theory  of  relativity  are  invoked.   Then, 
the  rate  of  change  of  both  mechanical  momentum  of  the  matter 
within  S  and  electromagnetic  momentum  of  the  field  within  V, 
may  be  equated  to  the  right-hand  side  of  eq.  (16)  which  must 
now  be  interpreted  as  the  inward  flow  of  momentum  through 
the  surface  S  per  unit  of  time  [Ref.  3] 


§£  (Pm  +  Pe)   =   /   [e(E.n)E  +  u(H-n)H  -  \    (eE2  +  yH2)n]  dS 

S  (17) 
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Transforming  this  surface  integral  into  a  volume 
integral  and  using  Maxwell's  equations,  the  resultant 
expression  would  have  the  form 

|.C?     +P    )    =    HpE  +  J  xB  -^(E2V£  +H2Vy)   +|T(DxB)]dv 
at     m        e  .7  2  dt 

V  (18) 

This  expression  gives  the  net  total  force  acting  on 
the  volume  V;  it  is  equal  to  the  rate  of  change  of  the 
total  momentum,  electromagnetic  and  mechanical. 

The  forces  exerted  by  the  field  on  charges  and 

neutral  matter  give  rise  to  an  increase  in  the  mechanical 

momentum  so  that  F  =  d/dt  P  .   This  force  can  be  determined 

m 

if  the  right-hand  side  of  Eq.  (18)  can  be  split  into  two 

parts  to  be  identified  with  P   and  P  .   However,  the  analysis 
c  m      e  2 

of  the  problem  may  become  difficult  since  this  separation 
is  not  an  obvious  task  to  carry  out.   That  is  why  the 
literature  on  the  subject  shows  that  several  hypotheses  have 
been  suggested  to  define  what  part  of  the  total  momentum 
is  associated  with  the  strictive  forces  in  the  medium  pro- 
ducing mechanical  deformation. 

In  order  to  simplify  the  analysis  of  the  mechanical 
forces  arising  from  the  interaction  of  a  travelling  electo- 
magnetic  wave  with  a  material  medium,  it  will  be  assumed 
that  the  medium  is  homogeneous,  devoid  of  free-charges 
and  of  mechanical  elasticity,  and  that  steady-state  condi- 
tions apply. 
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Under  these  assumptions,  Eq.  (18)  reduces  to  an 
expression  for  the  rate  of  change  of  electromagnetic 
momentum  giving  two  force  components  which  are  represented 
by  the  terms  JxB  and  3/3t(DxB)  .   As  JxB  is  independent 
of  time  variations,  it  may  be  interpreted  as  the  force 
that  would  arise  in  an  equivalent  DC  field;  on  the  other 
hand,  3/9t(DxB)  can  be  interpreted  as  the  force  associated 
with  the  travelling  electromagnetic  wave. 

It  is  then  clear  that  to  calculate  the  mechanical 

forces  on  the  material  medium,  the  time  rate  of  change  of 

the  electromagnetic  momentum  has  to  be  known. 

5.   Forces  Due  to  Electromagnetic  Energy  and  Velocity 
Gradients  in  a  Material  Medium 

The  total  energy  density  due  to  the  electric  and 

magnetic  fields  is  given  by 


W  =  jUE2  +  yH2)  (19) 


This  expression  represents  the  instantaneous  energy  stored 
in  the  electromagnetic  field  per  unit  area.   According  to 
the  theory  of  relativity,  mass  m  and  energy  W  are  equiva- 
lent and  must  obey  the  relation 


m  -  -|  (20) 

c 


Also,  if  mass  m  moving  with  velocity  v  has  a  rest  mass 
m  ,  the  following  relation  must  be  valid: 
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m 
«»  =      r   -  (21) 

y/l    ~  (V/C)2 

The  binomial  expansion  of  the  right-hand  side  of  Eq.  (21) 
will  result  in  the  expression 


m  =  mr+IV  =  mr  +  mr  £(|)  2  +  £<£)  4  +  *§(§)  6  +  . . .  ]      (22) 


m,  ,  the  kinetic  mass,  is  a  function  of  the  velocity  v; 
it  vanishes  when  v  =  0. 

Assuming  that  Eq.  (20)  can  be  applied  to  moving 
electromagnetic  energy  and  its  equivalent  mass,  then  the 
energy  W  can  be  split  into  a  rest  energy  W   and  an  energy 

of  relative  motion  W, 

k 


W  =  W  +  W,  =  (m  +m.  )c2  (23) 

r   k     r   k 


Consider  now  a  uniform  plane  wave  with  an  angular 
frequency  oo  and  a  phase  shift  constant  S  in  a  medium  with 
conductivity  a.   The  phase  velocity  v  of  the  travelling 
wave  is  given  by 


1/2 
I 

\iz  Ml  +  (a/cue)  *  +  1 


v  =  £  =  ( Z  )  (24) 


3       / ,   2 

[VI  +  (a' 


12     2 
If  the  wave  comprises  an  energy  density  W  =  w-(eE  +  uH  )  , 

a   z 

then  the  rate  of  flow  of  energy  per  unit  area  is: 
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s  =  4  (eE2  +  yH2)v  =  W  v  =  EH  (25) 


The  magnitude  of  the  electromagnetic  momentum 

-*■-*■-*■  •+■■*■ 

density  p  =  DxB   =   ye(ExH)    associated  with   the  moving 


energy  W  is 
a 

W  v 

EH     EH     a 
p  =  yeEH  =  •*-> —  =  — *•  =  — «■ 
r   K      l/y£    2     2 

'        V       V 

o     o 

(26) 


where  v  =  l//ye. 

The  mass-energy  equivalence  principle  applied  to 

2 

the  energy  W  would  be  expressed  by  the  relation  m  =  W  /v  . 
a  a    a  o 

Note  that  the  energies  W   and  W,  must  not  be  considered 

a      k 

the  same  quantities.   W,  is  an  additional  energy  of  motion 
for  a  rest  energy  W  and  it  translates  to  equivalent  mass 
in  terms  of  c  . 

For  a  transverse  electromagnetic  wave  progressing 
in  the  +z  direction,  the  electromagnetic  momentum  associated 
with  an  element  of  thickness  <5z  and  unit  cross-sectional 
area  is  given  by  the  expression 


W  v 

3p  =  -|-  6Z  (27) 

v 
o 


As  force  is  the  time  rate  of  change  of  momentum,  the 
corresponding  elementary  force  would  be 


A  1      dW         A 

a   -  &(«P>  =  4(vaE&  +  Kag)««  .(28) 

V 

o 
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From  basic  calculus 


d    9_  3_Z  _    d_ 
dt  ""  3Z  3t  '  V  3Z 


so  that 

6F  =  (Z-  -JL  +   a  |v}  62  ( 

v2  8Z    v2   9Z 
o        o 

The  integration  of  Eq.  (29)  allows  then  the  computation 
of  the  total  force  per  unit  surface  area  and  if  the  time 
averages  of  the  force,  F,  and  energy,  W,  are  considered, 
the  expression  for  the  total  time  average  force  per  unit 
surface  area  will  be 

2     aw  ,  a 

F  =  h  fjT  dZ  +  TT  /57av  If  dZ  (30) 

Eq.  (30)  denotes  the  presence  of  energy  and  velocity 
gradients  and  since  mechanical  forces  can  only  operate  on  a 
material  medium,  these  energy  gradients  must  arise  either 
from  dissipation  producing  heat  or  from  a  progressive 
space  change  in  the  stored  energy.   Also,  a  space  change 
in  the  propagation  velocity  yields  a  force  component  as  it 
is  the  case  when  power  is  reflected  at  an  interface  between 
two  different  media. 

As  stated  before,  the  gradients  of  energy  in  a 
source-free  medium  resulting  in  a  rate  of  change  of  the 
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electromagnetic  momentum,  may  be  produced  in  two  distinct 
ways: 

(i)   by  conversion  of  energy  to  heat, 
(ii)   by  transfer  of  stored  energy  to  or  from  the 
material  medium. 
In  the  source-free  conditions  assumed,  any  additional 
energy  put  into  the  material  as  the  wave  propagates  through 
it,  whether  to  provide  for  heat  losses  or  a  change  in  the 
stored  energy,  must  come  from  the  empty  space  which  the 
material  occupies  since  conservation  of  energy  requires  a 
constant  value  of  its  total  content.   Any  such  transfer  of 
energy  from  space  to  material  medium  is  associated  with  a 
force  acting  in  the  +Z  direction,  that  is,  it  is  a  positive 
force.   Then,  the  dissipation  of  energy  as  heat  or  a 
progressive  increase  of  the  stored  energy  must  produce  a 
positive  force. 

The  resultant  force  produced  in  a  material  medium 
through  the  process  described  is  called  radiation  pressure. 
6.   Energy  Gradients  —  A  More  Detailed  Analysis 

For  the  conditions  assumed,  it  follows  from 
Maxwell ' s  equations  that 


8E       9B 
x 


3z     at 


(3i; 


9H  dD. 

az*  "  »  Ex  +  at5  (32) 
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12      2 
Also,  the  expression  W  =  -^(eE   +  uH  )  may  be  written 

3     »-     X        V 


as 


Wa  =  I(£Ex  +  yHy)  (33) 


Combining  Eqs.  (31)  and  (32)  with  the  partial  derivative 
of  Eq.  (33)  with  respect  to  Z,  an  expression  for  the 
gradient  of  the  total  energy  in  the  electromagnetic  field 
can  be  readily  obtained 


9B  3D    i       o    n  o    ^ 

=  -uaE  H   -  £E   ^  -  yH    —2.  +  ±(E^  4t  +  H   4#) 
M   x  y      x  9t      y   3t    2V  x  3Z    y  3Z 

(34) 


Under  the  assumptions  being  considered  here,  i.e.,  charge- 
free  medium  and  unit  cross-sectional  area,  it  should  be 
noted  that  the  right-hand  side  of  Eq.  (34)  corresponds  to 
Eq.  (18)  and  therefore  it  gives  the  rate  of  change  of  the 
total  momentum  including  the  part  arising  from  strictive 
forces  producing  mechanical  deformation  which,  however,  is 
assumed  to  be  negligible  in  the  present  discussion. 

Since  the  medium  is  homogeneous  the  last  term  on 
the  right  side  of  Eq.  (34)  vanishes. 

When  the  fields  have  a  sinusoidal  time  variation, 
which  is  normally  the  case,  the  time  average  of  the  energy 
gradient  corresponding  to  Eq.  (34)  is  easily  computed  using 


31 


complex  wave  functions.   Thus,  representing  the  field 
components  as  complex  amplitudes  including  both  the  phase 
factor  and  the  space  dependence  of  the  wave  function,  the 
result  is 


3W  ,  ,  9H  8E 

__al  =  -yaRe(^ExH*)  -  y£Re(E*  -g.)   -  ueRe(H*  _£)      (35) 


where  Re  and  the  asterisk  stand,  respectively,  for  the 
real  part  and  complex  conjugate  of  the  quantities.   But 
noting  that 


3H  8E 

ReCEx  3t*>    =  "  Re(Hy  It**  (36) 


it  happens  that  Eq.  (35)  reduces  to  the  form 


8Wal      „  ,1 


dZ 


-yaRe(±-  E  H*)  (37) 

K     2   x   y 


This  time  average  of  the  gradient  of  the  total  energy  per 
unit  volume  is  associated  with  conversion  to  heat,  a  fact 
that  is  a  consequence  of  the  finite  conductivity  a   of  the 
medium;  as  this  gradient  is  a  negative  quantity,  it  repre- 
sents a  progressive  transfer  of  energy  to  the  material  with 
a  gradient  +yyaRe (E   H*)  and  so  the  corresponding  force  is 
positive. 

Eq.  (37)  does  not  reveal  any  gradient  of  stored 
energy  in  the  material;  therefore  it  is  necessary  for 
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further  analysis  to  show  the  possibility  of  the  existence 
of  such  a  gradient  in  the  medium  responsible  for  a 
supplementary  force. 

Representing  the  permittivity  and  permeability  of 
free  space  respectively  by  e   and  u  ,  the  instantaneous 
energy  stored  in  the  material  medium  alone  per  unit 
volume  is  given  by 


Wa2  "  I(£"eo)Ex  +  J'n-Vy  (38) 


Differentiating  Eq.  (38)  with  respect  to  Z,  taking  its 
time  average  using  complex  notation  and  assuming  a  homo- 
geneous medium  with  \i  —   y  ,  the  final  expression  for  the 
gradient  of  stored  energy  would  be 

8W  ,  ,     BE 

JzT-=   y(e-,o)Re(^H*  ^)  (39) 

In  the  case  of  a  pure  travelling  wave  this  gradient  will 
be  negative  [Ref.  5]  because  as  the  field  progresses 
through  the  medium  it  attenuates  due  to  its  finite  conduc- 
tivity.  The  gradient  represents  then  energy  progressively 
coming  out  of  the  medium  and  the  force  produced  will  be 
negative. 

So,  in  general,  it  appears  that  a  force  component 
may  result  from  a  rate  of  change  of  electromagnetic 
momentum  associated  with  stored  energy. 
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7.   Calculation  of  Radiation  Pressure  in  an  Infinite 
Homogeneous  Medium 

Eq.  (30)  is  the  basic  equation  for  the  calculation 
of  the  forces  produced  in  a  medium  with  characteristics 
y,  z   and  a,    due  to  a  rate  of  change  of  electromagnetic 
momentum  associated  with  a  travelling  wave  moving  in  the 
positive  Z  direction  with  a  constant  velocity  v.   Since 
the  velocity  is  constant,  the  second  term  on  the  right- 
hand  side  of  Eq.  (30)  is  zero  and  so  this  equation  becomes 
for  a  thickness  d, 


v2   d  8Wa 

*  "  T"  f      aT  dz  (40) 

v   0    3Z 
o 

3W  /9Z  can  be  decomposed  into  two  gradients  of  energy 
a 

/ 

according  to  Eqs.  (37)  and  (38)  so  that,  using  the  appro- 
priate signs  to  yield  positive  forces,  3W  /3Z  can  be  written 

a 

as 


8Wa     3W  ,    aw 

— i  =  -  — S±  +  — 2£  (41) 

SZ       3Z      3Z  ; 


and  Eq.  (40)  becomes 


F  =  F,  +  F2 


2d,  ,    3E 

=  ^j  f    [yaRe(|ExH*)  +y(£-£o)Re(|H*  ^)]dZ     (42) 

o 
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where  F.  and  F_  are  the  average  forces  resulting  from  the 
gradients  of  energy  associated,  respectively,  with  dissi- 
pation as  heat  and  with  stored  energy. 

a.   Forces  in  a  Medium  Extending  from  Z=0  to  Z=°° 

For  a  finite  conductivity  medium  extending  from 
Z=0  to  Z=°°,  the  total  power  associated  with  a  wave  is 
absorbed  which  means  that  the  overall  stored  energy  is 
zero.   Thus,  no  net  component  of  force  F2  can  be  expected. 
When  this  situation  occurs,  Eq.  (42)  takes  the  form 

2  oo 

F   =   F1   =  ±2     f      uaRe(j  Ex  H*)    dZ  (43) 

v      0  y 

o 

The  quantity  Rety  E   H*)  is  the  time  average  value  of  the 
instantaneous  power  density  s   carried  by  the  wave  which 
attenuates  with  Z  at  a  rate  determined  by  a,  the  attenuation 
coefficient  of  the  medium;  it  can  be  shown  [Ref.  6]  that 
s  is  given  by 

E2 

s   =  Re(i  Ev  H*)  =  ^t^T  cos  9  e~2ctZ  =  s,e~2aZ  (44) 
v2   x   y'    2  |r|  I  1 

where  E    is  the  value  of  E   at  Z  =  0  and  I n I  is  the  magni- 
xo  x  '  ' 

tude  of  the  complex  intrinsic  impedance  of  the  medium 

n  =  | n |  /e. 

Then  the  force  coming  from  energy  dissipation 
is  given  by 


2        °°  2         s 

F  =  Ft    =  2y     /      ya?1e"2aZ   =  ^vo  (^)  (45) 

o  o 
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and   since   v  =   2cc/ya 

^   =  4^  (46) 

v 
o 

b.   Forces  Due  to  Reflection  in  the  Medium 

If  the  wave  being  propagated  through  the  medium 
finds  a  partially  reflecting  surface,  then  a  certain  per- 
centage s  =  Wv  of  the  power  it  carries  and  its  associated 
electromagnetic  momentum,  suffer  a  complete  reversal  and 
so  the  second  term  on  the  right  side  of  Eq.  (30)  is  no 
longer  zero.   This  means  that  a  force  of  reflection  will 
be  produced 


v   -v       °  V 

o  o 


For  the  particular  case  of  a  lossless  medium  (a  =  0)  , 

the  fact  that  v  =  v   allows  a  further  simplification  of 

o  r 

Eq.  (47)  and  since  y  =  y   y   and  e  =  e   e 


o   r 


2s    2s 


F3  "  T-  '   ¥•  /£r  «r  (48) 

o 


c  is  the  velocity  of  light  in  free  space,  e   and  y   are 
respectively  the  dielectric  constant  and  relative  permea- 
bility of  the  medium  outside  the  reflecting  surface. 

As  seen  from  Eq.  (48) ,  the  force  due  to  reflec- 
tion is  directly  proportional  to  /e~~  ,  the  index  of  refrac- 
tion of  the  medium. 
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8 .   Rate  of  Change  of  Electromagnetic  Momentum  for 
A  DC  Field  and  its  Corresponding  Force 

When  frequency  approaches  zero,  the  velocity  v 
approaches  zero  too  and  the  same  happens  to  the  electro- 
magnetic momentum  as  can  be  seen  from  Eq.  (26) ,  for  example 

Consider  the  material  medium  immersed  in  a  DC 
field.   If  a  rest  energy  W,  analogous  to  W   is  considered, 
it  may  be  translated  to  equivalent  mass  m,  using  the 
velocity  c. 

As  a  logical  extension  to  this  case  of  'DC  flow 
of  energy ' ,  momentum  may  be  then  expressed  by  the  relation 


W     W 

p  «  nLC  -(-§)-  -£  (49) 

c 


So,  even  in  the  case  of  a  DC  field,  the  concept  of  an 
electromagnetic  momentum  responsible,  when  absorbed  by  the 
medium,  for  a  mechanical  force,  is  still  valid.   Without 
going  in  much  detail,  it  seems  natural  to  accept  the  idea 
of  a  smooth  transition  from  wave  propagation  to  DC  trans- 
mission.  It  should  be  noted  that  this  electromagnetic 
momentum  density  concept  is  equivalent  to  the  one  described 
by  Eq.  (4) . 

It  will  be  shown  next  that  under  DC  conditions,  the 
production  of  force  on  the  material  medium  via  conversion 
of  electrical  energy  to  heat  is  still  maintained. 

W,  /c  represents  the  volume  density  of  electromagnetic 
momentum  for  a  DC  field  in  a  homogeneous  medium.   When  the 
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field  components  are  transverse  to  the  direction  Z,  the 
electromagnetic  momentum  associated  with  an  element  of 
thickness  6Z  and  unit  cross-sectional  area  is 


<5p  =  _£  6z  (50) 


so  that  the  corresponding  elementary  force  is 


W        dW 
«F  =  k(«p)  =  |r(-|  «Z)  =  g^  6Z  (51) 


By  integration,  the  total  force  per  unit  surface  area  is 
given  by 


3W 
F  =  /  — -  dZ  (52) 

3Z 


For  time- invariant  field  components,  the  displace- 
ment current  density  3D  /3t  vanishes  and  so  Eq.  (32) 
becomes 

3H 


3Z   =  ^  Ex  (53) 


Using  Eqs .  (31)  and  (53) ,  Eq.  (34)  reduces  to 


3W 

*r  "  -"0EA  (54) 
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For  a  thickness  d  of  the  medium,  the  force  per 
unit  area  arising  from  dissipation  of  power,  according  to 
Eq.  (52)  ,  is  given  by 


F  =  mo      f      EH  dZ  (55) 

n     *  Y 


Obviously,  no  gradient  of  stored  energy  contri- 
butes to  this  mechanism  of  production  of  force. 
9 .   Summary  of  Important  Results 

The  discussion  of  the  forces  arising  in  a  material 
medium  under  the  action  of  an  electromagnetic  field,  was 
carried  out  in  a  simplified  version  and  results  in  a  few 
important  relations . 

It  was  shown  that  the  process  responsible  for  the 
production  of  such  mechanical  forces,  is  a  rate  of  change 
of  electromagnetic  momentum  associated  with  moving  energy 
or  its  equivalent  mass. 

Three  different  mechanisms  may  be  concurrent  in  the 
production  of  the  resultant  force,  each  one  being,  however, 
a  consequence  of  distinct  causes  originating  a  change  in 
momentum.   These  three  major  mechanisms  described  before 
are  repeated  here  along  with  the  corresponding  equations. 

a.   Change  in  Momentum  Due  to  Energy  Gradients 

Energy  gradients  in  the  medium  may  be  produced 
in  two  distinct  ways  giving  rise  to  different  force 
components : 
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(i)    By  conversion  of  energy  to  heat 


-   v2   d      1 

F  =  -j   /   yaRe(iExH*)  dZ  (56) 

v   0  y 

o 


(ii)   By  transfer  of  stored  energy  to  or  from 
the  medium 


-   v2    d  1    3EX 

F  =  7  o   y(e"Eo)Re(??  W^]    dZ  (57) 

o 

b.   Change  in  Momentum  Due  to  Velocity  Gradients 

A  force  component  can  also  arise  if,  for  a 

given  power  (W  v) ,  there  is  a  change  in  propagation 
a 

velocity 


F  =  -V/Wav  If  dZ  <58> 

vo 


B.   TRANSIENT  SURPACE  HEATING 

1.   Simple  Thermoacoustic  Source 

Another  process  of  the  conversion  of  electromagnetic 
to  acoustic  energy  to  be  described  now,  is  the  so-called 
thermal  stress  generation  by  transient  surface  heating  of 
an  elastic  medium. 

This  is  a  known  process  since  the  end  of  the  19th 
century  but  only  recently,  due  to  improvements  in  high  power 
laser  technology,  it  was  possible  to  make  an  experimental 
verification  of  the  complex  theory  governing  this  mechanism. 
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The  solution  of  the  related  problems  require  a 
certain  number  of  assumptions  so  that  the  mathematical 
models  used  apply  only  to  the  specified  conditions. 

The  transient  heating  of  the  surface  of  a  body  may 
be  produced  in  several  ways,  namely,  by  electron  bombard- 
ment or  by  electromagnetic  energy  absorption.   The  latter 
process  is  the  one.  to  be  considered  here,  that  is,  the 
conversion  to  an  elastic  wave  of  that  portion  of  the 
incident  electromagnetic  energy  absorbed  by  the  medium. 

It  will  be  adopted  a  one-dimensional  treatment 
of  the  problem  and  assumed  an  uniform  heating  at  and  near 
the  surface  of  a  semi-infinite  absorbing  medium  as  well 
as  an  incident  electromagnetic  intensity  uniformly  distributed, 

Under  these  conditions ,  a  temperature  gradient 
normal  to  the  surface  will  be  produced.   This  gradient  is 
characterized  by  the  physical  properties  of  the  medium  and 
by  the  nature  of  the  incident  flux.   The  consequent  ther- 
mal expansion  gives  rise  to  a  production  of  strains  in  the 
medium  which,  in  turn,  leads  to  the  generation  of  stress 
waves  propagating  away  from  the  heated  surface. 

The  time  profile  of  the  stresses  and  the  amplitude 
of  the  waves  produced  for  a  given  absorbed  power  density, 
depend  upon  the  acoustic  boundary  conditions  at  the  illum- 
inated interface,  that  is,  upon  the  elastic  constraints 
applied  to  the  heated  surface.   These  boundary  conditions 
may,  in  general,  be  divided  into  two  types:   stress-free 
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and  constrained  by  contact  with  another  body.   The  stresses 
and  wave  amplitudes  obtained  in  the  stress-free  condition 
may  be  smaller  than  those  obtained  under  constrained 
heated  surface. 

2 .   Equation  of  Motion 

Consider  that  the  plane  x  =  0  in  a  cartesian 
coordinate  system  is  the  surface  of  the  semi-infinite 
elastic  medium  being  illuminated  with  uniformly  distributed 
electromagnetic  energy.   Part  of  this  energy  will  be  re- 
flected at  the  interface  but  the  remaining  will  penetrate 
the  medium.   If  a  and  I   represent,  respectively,  the 
attenuation  constant  of  the  medium  and  the  radiation  inten- 
sity at  the  surface,  then  the  radiation  intensity  I(x)  at 
a  distance  x  from  the  plane  x  =  0  is  given  by  the  familiar 
exponential  law 


ICx)  =1   e"aX  (59) 

o 


This  means  that  the  electromagnetic  wave  is  being  absorbed 
by  the  medium  as  it  progresses  through  it. 

This  absorption  of  energy  causes  a  temperature  rise 
at  and  near  the  surface  and,  as  a  result  of  this  effect, 
a  strain  z      is  produced 


=  9u(x,t)  m   a9(X/t)  (60) 

X       dX 
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where  u(x,t)  is  the  x  component  of  particle  displacement , 
a  is  the  coefficient  of  linear  thermal  expansion  and 
9(x,t)  is  the  temperature  rise  above  the  uniform  initial 
temperature.   Recalling  that  the  problem  is  being  treated 
in  one  dimension 


9u(x,t)    3u(x,t)  _  n  ,--. 

— 3y SZ °  (61) 


therefore  the  strains  £   and  e„  are  also  zero. 

y     z 

In  order  to  produce  the  strain  e   in  the  absence 
of  heating,  a  stress  a  would  have  been  necessary 


a  =  -  Ba0  (62) 

x 


where  B  is  the  bulk  modulus  of  elasticity. 

In  the  presence  of  both  heating  and  stress,  it 
can  be  shown  [Ref .  7]  that  the  corresponding  stress-strain 
relationship  for  the  case  of  negligible  lateral  inertia 
and  shear,  exhibited  by  liquid  media  for  example,  is  given 
by 


°v  =  Bev  "  Ba6  =  B  ~  -  Ba9  (63) 

XX  d  X 


The  resulting  equation  of  motion  for  the  semi- 
infinite  body  is  then  expressed  in  the  form 
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d     U  _  X 

Tt7    9x 


k  (B  i  -  Bae)  <64> 


where  p  is  the  density  of  the  medium.   Noting  that  the 
modulus  of  elasticity  and  the  velocity  v  of  the  compressional 
wave  propagation  in  the  medium,  are  related  by 


B  =  pv2  (65) 


eq.  (64)  becomes 


2      2 

1   3  u    8  u    .,  96  tcc. 

■T  771  -  TT  "  a  3x  (66) 

V    dt      dX 


This  is  the  nonhomogeneous  wave  equation  to  be  solved  for 
the  elastic  particle  displacement. 
3 .   Temperature  Distributions 

In  order  to  solve  eq.  (66),  the  temperature  distri- 
bution function  8(x,t)  has  to  be  known. 

For  the  case  of  an  electromagnetic  wave,  if  it  was 
assumed  that  the  whole  energy  was  absorbed  in  the  plane 
x  =  0,  relatively  simple  approximate  temperature  distributions 
could  be  obtained.   However,  absorption  of  energy  actually 
takes  place  in  a  region  of  finite  thickness  which  is  related 
to  the  skin  depth.   If  this  fact  is  taken  into  account,  the 
temperature  distributions  are  obtained  through  a  more 
difficult  analysis  but  they  will  be  accurate. 
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A  particular  temperature  distribution  is  examined 
here:   it  is  the  one  resulting  from  a  constant  input  electro- 
magnetic intensity  of  short  time  duration,  suffering  an 
exponential  attenuation  in  the  medium.   It  is  assumed  that 
most  of  the  input  energy  degenerates  to  heat  and  that  the 
process  is  adiabatic. 

If  radiation  of  intensity  I0  hits  the  interface  in 
the  time  range  t=0tot=Tfca  heat  source  in  the  region 
x  >  0  will  be  created  having  a  heat  production  rate  deter- 
mined by  eq.  (59) .   The  resultant  heat  flow  can  be  considered 
one-dimensional  as  long  as  the  penetration  depth  a    is  much 
less  than  the  diameter  d  of  the  surface  receiving  heat. 

A  solution  for  this  temperature  distribution,  where 
the  heat  flow  is  supposed  to  be  in  the  positive  x  direction, 
was  achieved  by  Carslaw  and  Jaeger  [Ref.  8]  but  the  final  form 
of  the  function  9(x,t)  is  a  rather  complicated  expression. 

Under  certain  conditions,  the  solution  just  mentioned 
may  be  simplified.   In  fact,  an  approximation  can  be  found 
[Ref.  7]  for  the  case  of  intense  heating  of  media  having  low 
thermal  conductivity,  such  as  liquids,  and  for  durations  of 
the  radiation  pulses  from  10  to  50  nsec.   This  approximation 
is  expressed  by 

-ax 
a  I  te 

9(3t't)  -   gpc  (67) 

where  J  and  C  are,  respectively,  the  mechanical  equivalent 
of  heat  and  the  specific  heat  of  the  medium. 
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This  type  of  heating  can  be  produced  using,  for 

example,  a  Q-switched  laser  with  output  in  the  range 

2 
50-100  MW/cm  .   The  temperature  rises  so  obtained,  although 

rapid,  do  not  exceed  tens  of  degrees  but  the  experienced 

temperature  gradients  would  be  extremely  large.   As  an 

example  [Ref .  7] ,  it  could  be  expected  in  water  colored  by 

2   -1 
a  dye  to  give  a  =  10  cm   ,  temperature  gradients  as  high 

3 

as  7.5x10   °C/cm  and  temperature  rates  of  change  of  the  order 

9 

of  2.5x10   °C/sec  at  the  surface  during  the  laser  pulse. 

If  nonmetalic  media  is  considered,  the  temperature 
decay  for  t  >  T  is  a  slowly  varying  function  of  time, 
becoming  significant  only  for  t  greater  than  milliseconds. 
Also,  the  elastic  transients  produced  through  this  process, 
propagate  away  from  the  heated  place  in  periods  of  time  as 
short  as  microseconds.   So,  it  can  be  assumed  that  for  t  >  T 

al  Te_ax 
9(x't)  =-^c (68) 


Equations  (67)  and  (6  8)  represent  then  the  source 

functions  to  be  used  in  eq.  (66)  in  order  to  determine  the 

equations  for  displacement  and  stress. 

4.   Solution  of  the  Equation  of  Motion 
For  Two  Boundary  Conditions 

a.   Constrained  surface  (rigid  boundary) 

Substituting  eqs.  (67)  and  (68)  in  eq.  (66) , 

the  resultant  wave  equations  during  and  after  the  heating 

pulse  are 
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2         2 

3  u  =  v2  JTu  + 

at     ax2 


2T7.  -ax 
v  Wte      , 


2TTm  -ax 

I  v  WTe      , 


t  <  T  , 


t  >  T 


(69) 


where 


W  = 


2X 

aa  I 
c 

JpC 


(70) 


The  solution  of  these  equations  for  the  initial 
conditions 


3t 


=  u(x,0)  =  0 


(71) 


x,0 


and  the  appropriate  boundary  condition  for  a  constrained 
surface 


u(0,t)  =  0 


(72) 


would  give  the  expressions  for  the  particle  displacement 
u(x,t)  in  the  time  ranges  T  <  t  <  x/v  ,   x/v  <  t  <  x/v  +  T 
and  t  >  x/v  +  T   [Ref .  7] .   Once  u(x,t)  is  known,  3u/3x 
can  be  substituted  in  eq.  (63)  to  obtain  the  expression 
for  the  stress  in  the  medium  as  a  function  of  time  and 
distance  from  the  surface.   Since  the  condition  a    <<  d 
is  met,  an  elastic  plane  wave  is  then  launched  from  the 
surface  and  propagates  away  from  it. 
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Again,  due  to  the  extreme  complexity  of  the  equations 
for  displacement  and  stress  actually  obtained,  it  is  desirable 
to  try  some  sort  of  valid  simplification  to  render  them  more 
suitable  for  specific  applications.   This  is  possible  if  the 
values  of  x  are  chosen  such  that  ax  >>  1,  that  is,  the  points, 
x  should  be  those  far  from  the  source.   Then,  introducing 
a  reduced  time  variable  x  =  t  -  x/v,  the  approximate 
expressions  for  the  stress  would  take  the  form 


val 
cr(x't)  =  2JC2  " 


,,     -avTv  avx  ^rt 

(1  -  e    )e  ,   x<0  , 


[2  .  e-<*vT  _  eav(x-T)]    0<x<T^    (?3) 


,  avT    , x  -avx  __ 

^  ■  (e    -  1)  e  ,   x>T   . 


From  these  equations  it  can  be  observed  that  the 
maximum  stress  occurs  at  x  =  T/2  and  its  value  is  given  by 


,     =  !!J°  (1  -  e"-^2)  (74) 

max   2JC 


which  approaches  the  limiting  value  of  val  /2JC  for  avT/2  >>  1 
b.   Free  surface  (pressure-release  boundary) 

The  boundary  condition  at  the  stress-free  surface 
is  now 


aC0,t)  =  B  3u(^t}  -  Ba6(0,t)  =  0  (75) 

d  X 


or,  equivalently, 
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^§#  =  ae(0,t)  = 


aaIot 
JpC 


aal  T 
o 

JpC 


t<T 


t>T 


(76) 


Under  this  condition  and  following  a  procedure  similar  to 
the  one  used  for  the  constrained  surface,  it  is  possible  to 
get  a  set  of  approximated  equations  describing  the  stress  at 
observation  points  far  from  the  heated  surface,  ax  >>  1, 
and  in  terms  of  T 


val 


crCx,t)  = 


2JC 


,,     -avT.   avT 
f  (1  -  e     )  e 


r  -avx     av(T-T) , 
•  <  [e     -  e       ] 


,,     avT.  -avT 
L  (1  -  e    )e 


x<0    , 
,   0<T<T   , 

T>T 


(77) 


These  equations  tell  that  maximum  stresses  having  the  value 


val 


max   2JC 


o  ,,     -avT, 
(1  -  e     ) 


(78) 


appear  at  x  =  0  and  x  =  T  and  approach  a  limiting  value  of 

val  /2JC  for  avT  >>  1  . 
o 

5 .   Stress  Waves  in  the  Medium 

Solving  eqs .  (73)  and  (77)  for  various  values  of 
the  parameter  avT,  a  plot  of  stress  versus  reduced  time  can 
be  obtained. 
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Typical  curves  obtained  in  this  manner  would  look 
like  those  shown  in  Fig.  1. 


time  t 


(a) 


(b) 


Fig.  1.   Typical  stress-time  plots  for  constrained 
surface  (a)  and  stress-free  surface  (b) 


The  constrained  surface  situation  gives  rise  to  stress 
consisting  of  a  unidirectional  compressive  pulse,  while  in 
the  free  surface  case  the  stress  suffers  an  abrupt  change 
from  compression  to  tension  during  the  time  T. 

To  give  an  idea  about  the  theoretical  magnitude  of 
the  resultant  stresses,  the  use  of  the  free-surface  expressions 
for  the  case  of  a  100  MW,  30  nsec  laser  pulse,  illuminating 
water  with  a  =  100  cm   ,  would  lead  to  produced  stresses  with 
peaks  of  ±20  atm. 

As  expected  from  the  theory,  the  peak  amplitude  of 
the  pulse  decreases  and  its  duration  increases  as  the 
attenuation  constant  of  the  medium  is  reduced. 
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High  power  laser  pulses  may  be  then  responsible 
for  the  production  of  thermal  effects  in  media.   The  rapid 
changes  in  temperature  occuring  at  a  surface  absorbing 
rapidly  changing  amounts  of  energy  is  due  to  the  fact  that 
thermal  conduction  is  a  slow  process.   The  elastic  waves 
detected  at  points  far  from  the  surface  of  a  medium  as  a 
result  of  the  absorption  of  an  electromagnetic  pulse  of 
high  frequency,  will  have  a  frequency  content  consisting 
primarily  of  frequency  components  associated  with  the  pulse 
envelope  rather  than  with  the  carrier. 

A  few  experiments  on  the  generation  of  acoustic 
signals  in  liquids  by  laser  pulses,  had  been  performed 
recently  and  the  results  seem  to  confirm  the  theory  behind 
the  transient  heating  process  of  elastic  wave  generation. 

E.F.  Carome,  N.A.  Clark  and  C.E.  Moeller  [Ref.  9], 
using  a  laser  beam  from  a  0.1  Joules,  50  nsec  single  pulsed 
ruby  laser,  verified  the  validity  of  eqs.  (73)  and  (77)  for 
the  case  of  transient  heating  of  water  surface.   The  measure- 
ments were  performed  with  various  concentrations  of  Prussian 
blue  in  order  to  change  the  attenuation  constant  of  the  water 
The  shape  of  the  pulses  obtained,  closely  approximate  those 
predicted  by  theory.   As  observed,  an  increase  in  a,  that  is, 
an  increase  in  the  absorptivity  of  the  medium,  increased  the 
elastic  wave  amplitude.   This  result  indicates  that  the 
produced  elastic  waves  could  not  have  been  generated  through 
a  process  of  radiation  pressure  by  reflection  at  the  surface 
since  an  increase  in  absorptivity  decreases  the  amount  of 
reflected  energy. 
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6.   Efficiency  of  Energy  Conversion 

The  conversion  efficiency  N  of  this  process  is 
defined  as  the  ratio  of  the  energy  of  the  propagated  elastic 
wave  E   to  the  energy  of  the  incident  electromagnetic  wave  E 

E 
N  =  =4  (79) 

Ee 

As  may  be  guessed,  this  process  is  characterized  by 

a  very  low  efficiency. 

For  both  free  and  constrained  surface  conditions,  the 

energy  E   is  a  function  of  the  area  through  which  the  elastic 
a 

wave  is  propagated  and  the  area  over  which  the  electromagnetic 
wave  is  incident.   Assuming  that  both  areas  have  the  same 
value,  the  resultant  expressions  for  N  [Ref.  7]  are,  for  the 
case  of  free-surface 

va  I 

N   =   5-0-  •  F,  (avT)  C80) 

2p(TJ       L 


where 


F.CavT)   =  L^—^E* (81, 


and  for  the  case  of  constrained  surface 


va2I 


N   =  *-%   •  F9(ctvT)  (82) 

2Pcrcr   Z 
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where 


F2(avT)   =  H2Te +  3*avT  +  2otT  '  3       (83) 


In  these  expressions  I   is  in  watts  per  square  centimeter 
while  all  the  other  quantities  are  in  the  cgs  system  of  units. 

A  plot  of  the  functions  F,  and  F-  versus  avT,  allow 
an  easier  and  quicker  numerical  calculation  of  N.   These 
plots  show  that  the  maximum  value  of  F,  is  0.3  for  avT  =  2.0 
and  that  F~  approaches  a  maximum  of  2.0  for  large  values  of 
avT. 

As  an  illustration  of  a  numerical  calculation  of  the 

efficiency  showing  its  dependence  on  the  properaties  of  the 

medium,  the  value  N  =  1.83x10    I  F..  (avT)  is  obtained  for 

o  1 

-9 
water  and  N  =  6.14x10   I  F, (avT)  for  carbon  tetrachloride, 

o  1 

both  for  the  free  surface  condition. 

A  final  remark  about  the  functional  dependence  of  N 

on  the  variable  I   should  be  made.   From  the  anaysis  of 

eqs .  (80)  and  (82)  it  seems  apparent  that  efficiency  can  be 

increased  without  limit  as  I   is  increased.   However,  this 

o 

cannot  be  a  true  statement  since  beyond  a  certain  value  of 
input  power  the  assumed  model  of  stress  generation  is  no 
longer  valid.   Some  limitations  of  this  model  result  from 
the  following  facts.   First,  the  thermal  characteristics  of 
the  medium  impose  limits  on  the  radiation  intensity  since 
phenomena  like  boiling  in  liquids  or  melting  and  vaporization 
in  solids  may  occur.   Stresses  still  continue  to  be  generated 
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although  by  a  different  process.   Secondly,  since  at  high 
intensities  the  thermal  and  elastic  properties  of  the  medium 
become  temperature  dependent,  the  formulas  used  before  are 
no  longer  applicable.   And  finally,  for  significant  values 
of  efficiency,  an  appreciable  fraction  of  the  input  energy 
is  lost  by  radiation  and  so  the  temperature  distribution 
assumed  cannot  be  used. 

7.   Thermoacoustic  Arrays 

The  way  the  conversion  mechanism  has  been  treated 
here  through  a  simplified  thermodynamic  model  emphasizing 
the  process  of  conversion  of  laser  light  to  acoustic  energy 
in  liquids,  may  be  considered  an  engineering  approach  of 
great  practical  interest. 

This  phenomenon  can  be  analyzed  in  a  more  general 
sense  to  include  the  treatment  of  a  wider  range  of  specific 
cases  with  different  types  of  media  and  electromagnetic 
sources . 

R.M.  White  [Ref.  10]  in  a  theoretical  and  experimental 
work  analyzed  the  conversion  process  for  several  temperature 
distributions  due  to  different  input  energy  fluxes. 

The  generation  of  low  frequency  sound  from  the 
absorption  of  light  has  become  a  subject  of  growing  importance 
with  potential  significance  in  underwater  acoustics. 
Westervelt  and  Larson  [Ref.  11]  in  a  recent  theoretical 
formulation,  predicted  the  development  of  highly  directive 
sound  beams  in  water  by  means  of  thermalization  of  modulated 
laser  light.   The  generation  of  this  highly  directional 
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acoustic  wave  results  from  an  extension  of  the  concept  of 
a  simple  thermoacoustic  source  to  that  of  a  thermoacoustic 
array,  that  is,  a  series  of  thermoacoustic  sources  where 
pressure  fields  combine  to  give  a  directional  effect. 

A  theoretical  summary  of  the  mechanism  leading  to 
the  preduction  of  such  a  wave,  which  propagates  perpendicular 
to  the  axis  of  the  laser  beam,  is  given  next. 

Assume  that  a  laser  radiates  an  amplitude  modulated 
beam  of  cross  sectional  area  S   propagating  in  the  +z 
direction  as  shown  in  Fig.  2. 


P(r,0) 


Fig.  2.   Geometry  of  Laser  Excitation 

Due  to  the  absorption  of  energy  in  the  medium,  each 
elemental  volume  experiences  a  rise  in  temperature  which, 
in  turn,  causes  an  oscillatory  expansion  in  the  thermal 
volume.   If  the  beam  illuminates  the  surface  of  a  medium  of 
low  thermal  conductivity,  the  solution  of  the  wave  equation 
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for  the  pressure  distribution  resulting  from  the  exponential 
absorption  of  energy  is 


jaa)aIo      -jut         r      -az        eja)/v'r-r'i     _3    ,         ,Q., 
P      =      -      47TC        e  /e  — -^— d   r»  (84) 

r-r ' 


where  co  is  the  laser  angular  modulation  frequency,  a  the 
logarithmic  coefficient  of  thermal  expansion,  C  the  specific 
heat  per  unit  area,  I   the  output  intensity  of  the  laser, 
a  the  attenuation  coefficient  for  the  laser  beam  and  v  the 
sound  velocity.   This  result  was  obtained  for  a  source 
function  representing  the  energy  per  unit  volume  per  unit 
time  absorbed,  given  by   I  exp(-az-jwt) . 

Looking  for  a  far  field  approximation,  r  has  to  obey 
the  condition 


r»-«-  =  4  (85) 

a  v     a 


So,  simplifying  and  performing  the  integraion  on  eq.  (84), 
the  acoustic  pressure  field  radiated  by  this  mechanism 
turns  out  to  be 

jaa)aPQ    ejkr-jajt 
p  =     4-rrrC   *  a  +  jkcos  9  (86) 

where  P   is  the  laser  output  power.   Averaging  eq.  (86) 
over  one  period,  the  resulting  acoustic  intensity  would  be 
given  by 
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,  ao)aPQ   2  , 

1      =      2pT    *     (47rrC    ]       *       2^,2         27  (87) 

K  a      +  k   cos    0 


In  this  equation  p  represents  the  density  of  the  medium. 

The  maximum  value  of  eq.  (87)  occurs  at  0  =  tt/2 
and  the  angle  at  which  the  intensity  is  reduced  to  one  half 
its  maximum  is 


'1/2   -   f  *  f  (88> 


The  total  acoustic  power  generated  when  the  limit  condition 
k/a  >>  1  applies,  is  given  by 

2   2 

acoa  P 

w  =   2_  (89) 

16pC 

The  radiation  pattern  given  by  eq.  (87)  is  symmetric  about 
the  Z  axis  and  has  a  sharp  maximum  at  0  =  tt/2  in  the  limit 
a  <<  k  . 

Equations  (88)  and  (89)  show  that  when  a  increases 
giving  rise  to  more  intense  thermal  expansions,  larger 
acoustic  power  outputs  are  produced,  but  when  a  decreases, 
the  length  of  the  array  increases  and  so  the  directivity 
is  improved. 

On  the  other  hand,  for  a  given  value  of  a,  increasing 
the  modulation  frequency  co,  the  beam  may  be  made  as  narrow 
as  desired.   An  increase  in  conversion  efficiency  is  also 
achieved  as  oo  is  increased. 


57 


The  array  acts  then  like  a  broadside  radiator 
where  directivity  is  proportional  to  the  attenuation 
coefficient. 

All  these  theoretical  results  were  verified  and 
validated  by  experiments  performed  by  T.G.  Muir,  C.R.  Cul- 
bertson  and  J.R.  Clynch  [Ref .  12] .   These  experiments  showed 
further  that  for  applications  of  laser  excited  arrays  to 
ocean  acoustics,  high-power  lasers  in  the  green-blue  region 
are  required  in  order  to  obtain  extremely  narrow  acoustic 
beams . 

C.   STIMULATED  BRILLOUIN  SCATTERING 
1.   Brillouin  Scattering 

Sound  waves  in  a  medium  may  also  be  created  from 
the  interference  of  two  intersecting  laser  beams.   This 
effect  has  been  verified  experimentally  [Refs.  13,14]  and 
its  detailed  theoretical  formulation  may  be  found  in  the 
literature  [Refs.  15,16]. 

A  scattering  phenomenon  can  be  defined  for  the 
purposes  of  the  present  discussion  as  a  process  whereby 
a  physical  system  extracts  energy  from  a  beam  of  light  and 
re-emits  this  energy  at  the  same  or  at  a  different  wave- 
length, in  the  same  or  different  directions  from  the  incident 
beam. 

Scattering  may  occur  in  a  number  of  different 
situations  but  among  all  the  forms  it  can  take,  Brillouin 
scattering  is  the  one  to  consider  as  far  as  conversion  of 
electromagnetic  to  acoustic  energy  is  concerned.   This  type 
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of  scattering  is  an  interaction  process  between  acoustic 
waves  at  microwave  frequencies  and  electromagnetic  waves 
at  optical  frequencies. 

It  has  been  accepted  that  an' electromagnetic  field 
has  both  a  particle  and  a  wave  nature.   Such  a  statement 
leads  to  the  definition  of  photon  as  a  particle  charac- 
tenzed  by  zero  mass,  energy  hv  and  wave  vector  k.   Simi- 
larly, the  treatment  of  lattice  modes  of  vibration  in 
material  media,  leads  to  the  phonon  concept,  a  particle 
having  zero  mass,  energy  hv  and  momentum  hg  where  v  is 
the  phonon  frequency,  g  the  phonon  wave  vector,  h  =  h/2ir 
and  h  the  Planck's  constant. 

The  acoustic  waves  involved  in  the  production  of 

6       12 
Brillouin  scattering,  usually  range  from  10   to  10   Hz 

and  can  be  externally  introduced  into  a  material  medium, 

can  be  already  present  in  the  material  in  the  form  of 

thermally  activated  phonons  or  can  be  stimulated  in  the 

medium  by  the  light  beam  itself. 

2.   Bragg  Reflection 

Consider  that  a  sound  wave  at  angular  frequency 

u)   is  traveling  upward  in  a  cubic  block  of  a  material 

medium,  as  shown  in  Fig.  3  in  a  cross  section  perspective. 

When  a  plane  wave  of  light  of  angular  frequency 

03   hits  the  left  face  of  the  block  in  a  direction  such 
o 

that  its  wavefronts  are  perpendicular  to  the  acoustic 
wavefronts,  an  interaction  process  will  take  place  within 
the  material. 
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FIG.  3  —  The  Debye-Sears  Effect 

The  results  of  this  scattering  process,  called 

Debye-Sears  effect,  is  similar  to  what  happens  when  a 

carrier  frequency  co   is  phase  modulated  with  a  single  tone 

co  /  that  is,  light  at  angular  sideband  frequencies  co  , 
m  o 

co  +co  #  co  ±2co  ,  etc.,  will  emerge  from  the  right  face  of  the 
o  m   o   m  3  3 

block.   If  X  and  A  represent,  respectively,  the  light  and 

sound  wavelengths,  then  the  wavefronts  associated  with  the 

first  upper  and  lower  sidebands  co  ±co  ,  will  be  tilted  by 
rr  o   m  * 

an  angle  9  z    ±X/A  with  respect  to  the  undeflected  carrier 
light  co  .   Similarly,  for  the  nth  order  sidebands,  the 
deflection  angle  will  be  ±  nX/A. 
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The  output  intensity  of  the  various  sidebands  is 

the  same  not  only  because  their  amplitudes  have  to  obey 

Bessel  functions  behaviour,  but  also  because  a  destructive 

interference  occurs  in  the  medium.   For  a  fixed  co  ,  the 

m 

extent  in  which  each  sideband  is  'destroyed'  depends 
upon  its  frequency  as  well  as  upon  the  physical  dimensions 
of  the  material.   So,  this  type  of  scattering  will  diffract 
just  a  small  amount  of  the  total  input  light. 

Suppose  now  that  the  input  ray  is  rotated  to  make 
an  angle  a  ~  X/2A  with  the  acoustic  wavefronts.   The  process 
developed  is  equivalent  to  that  occuring  when  the  incident 
ray  was  parallel  to  the  acoustic  wavefronts. 

As  the  first  upper  sideband  still  makes  an  angle 
A/A  wit;h  the  undeflected  output  light,  the  interaction  pro- 
cess for  this  sideband  may  be  though  of  as  a  simple  reflection 
of  light  from  the  acoustic  wavefronts.   This  effect  is 
depicted  in  Fig.  4. 
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—  The   Bragg  Effect  Geometry 
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Moreover,  this  sideband  is  the  only  one  suffering 
a  constructive  interference  in  the  medium  and  therefore 
it  is  the  only  one  with  significant  intensity. 

The  angle  of  incidence  a  such  that  a  =  A/2 A,  or 
more  precisely 


a  =  sin"1  ^.  (90) 


is  called  Bragg  angle.   Clearly,  a  Bragg  angle  can  be 
defined  for  each  output  frequency. 

Due  to  this  mechanism,  the  Bragg  effect  allows  the 
selection  of  a  particular  mode  of  output  light  by  a 
simple  control  of  the  direction  of  the  incident  light, 
that  is,  by  forcing  light  and  sound  to  meet  at  the  appro- 
priate angle. 

It  can  be  shown  that  at  exactly  the  Bragg  angle, 
the  scattering  attains  its  maximum. 
3 .   Production  of  Sound  Waves 

When  two  light  beams  with  frequencies  go,  and  lo_ 
and  wave  vectors  K,  and  K„  intersect  in  a  medium,  an 
induced  electrostrictive  pressure  wave  of  frequency 
|go,-go2|  and  wave  vector  K,-K?  may  be  produced.   If  the  angle 
of  intersection  of  the  light  beams  is  such  that  the  phase 
velocity  co, -go2/ |K. -K2  |  of  the  pressure  wave  is  the  same 
as  the  velocity  of  the  sound  in  the  medium,  then  a  traveling 
wave  will  be  produced.   It  is  verified  that  the  angle 
satisfying  this  condition  is  precisely  twice  the  Bragg 
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angle.   It  should  be  noted  that  this  effect  can  only  occur 
for  an  input  beam  energy  above  a  well  defined  threshold. 

This  mechanism  of  interaction  of  intense  light 
beams  leading  to  the  production  of  coherent  sound  is  called 
stimulated  -Brillouin  scattering. 

Consider  now  the  situation  shown  in  Fig.  4  where 
light  and  sound  intersect  at  the  Bragg  angle.   The  inci- 
dent and  scattered  light  wavefronts  intersect  at  twice  the 
Bragg  angle  and  so  they  should  form  an  interaction  pattern 
giving  rise  to  stimulated  Brillouin  scattering.   Thus, 

additional  sound  at  u  is  launched  from  the  interaction 

m 

region. 

It  turns  out  that  in  the  case  of  Fig.  4(a),  the 
produced  (additional)  and  original  sound  waves  are  in 
opposite  phases  causing  a  progressive  decrease  of  the 
sound  amplitude  (deamplif ication) .   On  the  other  hand,  if 
the  situation  shown  in  Fig.  4(b)  occurs,  the  produced  and 
original  sound  waves  are  in  phase  and  amplification  of 
sound  is  then  verified. 

The  stimulated  Brillouin  scattering  mechanism  may 
be  reduced  to  the  following  picture:   if  the  input  energy 
is  high  enough,  instability  conditions  are  created  within 
the  medium  and,  even  when  no  input  sound  is  present,  pro- 
duction and  buildup  of  sound  at  to  =  |  oj. -w- j  starting  from 
thermally  activated  phonons,  will  be  verified.   Or,  in  other 
words,  the  elctrostrictive  mixing  of  two  light  beams  gives 
rise  to  amplification  of  thermal  sound  in  the  medium. 
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Brillouin  scattering  is  a  classic  example  of 
parametric  interactions  in  traveling  wave  systems  where 
the  power  transfer  between  the  waves  may  be  calculated  on 
the  basis  of  either  the  Manley-Rowe  relations  or  simple 
photon-phonon  energy  conservation  arguments.   It  is  impor- 
tant to  realize  that  the  case  of  stimulated  Brillouin 
scattering  gives  an  example  where  the  principle  of  para- 
metric oscillations  can  perform  a  useful  function,  namely, 
the  output  energy  at  the  lower  frequency  is  in  the  form 
of  sound  as  contrasted  to  the  electromagnetic  energy  of 
the  pump. 

It  has  been  proved  experimentally  [Ref.  13]  that 
the  observed  sound  cannot  be  the  result  of  thermal  shock 
in  the  medium.  In  fact,  removing  one  of  the  light  beams 
without  affecting  the  total  input  light  power,  it  can  be 
observed  that  sound  is  no  longer  produced.  That  is,  the 
signal  observed  at  the  beat  frequency  of  the  two  light 
beams  is  indeed  optically  induced  sound. 
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III.   EXPERIMENTAL  PROCEDURE 

A.   APPARATUS  DESCRIPTION  AND  RESULTS 
1.   Radiation  Pressure  Experiment 

In  the  first  experiment  performed  it  was  intended 
to  observe  the  forces  of  radiation  pressure  developed  in 
fresh  water  by  a  plane  electromagnetic  wave. 

Due  to  the  mismatch  between  air  and  water,  part 
of  the  energy  is  reflected  back  at  the  interface  and  part 
of  it  is  transmitted  through  the  water  which  is  assumed 
here  to  be  devoid  of  mechanical  elasticity  and  of  free 
charges. 

To  perform  the  experiment,  a  portable  radar  model 
AN/PPS-6 (XE-6)  was  used  as  a  source  of  electromagnetic 
field.   Some  characteristics  of  interest  of  this  radar 
are: 

pulse  repetition  rate 2  KHz 

frequency 9.3  GHz 

pulse  length 0.5  ysec 

peak  power 100  W 

average  power 47  mW 

power  supply  requirements 13  V  at  1.5A 

The  power  supply  chosen  was  a  POWER  DESIGNS  INC.,  model 
3240.   Appropriate  pieces  of  coaxial  cable  and  rectangular 
waveguide  both  assumed  lossless,  were  used  to  connect  the 
radar  output  to  a  horn  antenna,  a  SPERRY  MICROLINE  model 
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56x1.   According  to  its  calibration  curves  an  absolute 
gain  of  15.8  db  can  be  expected  at  9.3  GHz.   The  antenna 
was  oriented  to  face  the  surface  of  the  water  with  the 
axis  nearly  perpendicular  to  it  so  that  a  transverse  uni- 
form plane  wave  could  be  assumed  at  the  illuminated  surface 

In  order  to  detect  the  pressure  waves  produced  in 
the  water,  an  hydrophone  type  LC-10  was  used.   This  is  a 
small  device  about  5.4  cm  long  having  a  circular  cross 
section  with  diameter  about  0.96  cm.   The  useful  steady 
state  bandwidth  of  this  acoustic  detector  extends  from 
0.1  Hz  to  120  KHz.   The  typical  free-field  voltage  sensi- 
tivity within  this  range  lies  between  -105  and  -110  db 
relative  to  1  V/ybar.   The  actual  hydrophone  used  was  the 
LC10  serial  #2122  exhibiting  a  flat  free  field  frequency 
response  between  1  KHz  and  9  KHz  at  -110  db. 

The  free  field  voltage  sensitivity,  or  microphone 
sensitivity  level  ML,  is  defined  by  the  relation 


ML  =  20  log  - .- /  c 

3  1V/ y bar 


where  V  is  the  amplitude  of  the  voltage  output,  P  is  the 
pressure  amplitude  of  the  acoustical  field  and  lV/ybar  is 
the  reference  sensitivity. 

The  LC10  typical  directivity  pattern  in  the  plane 
perpendicular  to  the  axis  is  almost  circular.  During  the 
measurement,  this  device  was  held  in  the  water  in  a  posi- 
tion parallel  to  the  surface  and  in  the  same  vertical  of 
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the  transmitting  horn  antenna.   The  hydrophone  depth  could 
be  varied  since  it  was  bonded  to  a  horizontal  plastic  rod 
45  cm  long  which  in  turn  was  attached  to  a  vertical 
metallic  rod;  this  assembly  could  be  moved  vertically  in 
front  of  a  tape-measure  stuck  on  a  fixed  wall. 

Through  a  25  ft.  coaxial  cable  the  hydrophone 
output  signal  was  fed  to  a  HEWLETT  PACKARD  voltage  ampli- 
fier model  465A  capable  of  either  20  or  40  db  of  amplification. 

Two  different  schemes  were  tried  to  measure  the 
forces  of  radiation  pressure.   In  the  first,  a  cubic  wood 
box  about  6  ft.  high  was  placed  over  one  of  the  acoustic 
tanks  in  Spanagel  Hall  Room  025.   The  box  was  covered 
inside  with  a  good  absorber  of  electromagnetic  energy;  ' . 
the  transmitting  equipment  was  placed  on  its  top  side 
where  a  small  hole  was  made  to  let  the  horn  antenna  go 
inside.   The  hydrophone  assembly  was  stuck  to  one  of  the 
vertical  wood  walls.   The  distance  between  the  horn  and  the 
surface  of  the  water  was  fixed  at  about  1  m. 

The  output  of  the  voltage  amplifier  was  connected 
to  the  input  of  a  TEKTRONIX  oscilloscope  type  422  either 
directly  or  through  a  very  narrow  band-pass  filter.   This 
filtering  action  was  achieved  through  the  operation  as  a 
narrow  band-pass  filter  of  a  wave  analyzer  HEWLETT  PACKARD 
model  302A. 

Due  to  the  considerable  environmental  random  noise 
amplitude  and  to  the  extreme  weakness  of  the  signal  pro- 
duced by  the  energy  conversion  mechanism,  no  effect  could 
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be  observed  in  the  oscilloscope.   No  signal  component  at 
about  2  KHz  appeared  when  the  output  signal  was  passed 
through  the  band-pass  filter,  neither  any  signal  resembling 
the  envelope  of  the  transmitted  waveform  could  be  seen 
over  the  noise  when  the  output  of  the  voltage  amplifier 
was  fed  directly  to  the  oscilloscope.   The  magnitude  of 
the  signal  to  be  expected  is  computed  later  in  the  theor- 
etical computations. section. 

In  the  second  scheme  used  it  was  intended  to  take 
advantage  of  a  crosscorrelation  or  autocorrelation  process. 

Instead  of  making  the  detection  in  the  big  acoustic 
tank  mentioned  before,  a  small  portable  water  container 
was  used.   It  was  covered  inside  with  styrofoam  to  minimize 
acoustic  reflections;  its  internal  dimensions  were  approxi- 
mately 26x30x67  cm.   This  time,  the  horn  antenna  was 
held  at  about  10  cm  above  the  surface  of  the  water. 

A  tunable  envelope  detector  was  installed  in  the 
wave-guide  in  order  to  get  the  envelope  of  the  transmitted 
waveform  to  be  crosscorrelated  with  the  received  signal  via 
hydrophone.   A  correlation  and  probability  analyser  SAICOR 
model  SAI-43A  with  capability  of  performing  either  auto  or 
crosscorrelation  seemed  to  be  a  suitable  equipment  to  use 
in  front  of  the  amplified  signal.   The  output  of  this 
analyzer  may  be  displayed  in  any  ordinary  oscilloscope. 

When  the  autocorrelation  of  the  envelope  of  the 
transmitted  waveform  was  tried,  no  result  was  obtained 
since  the  detected  pulse  was  so  narrow  that  even  in  the 
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fastest  sampling  mode  of  the  analyser  no  autocorrelation 
could  be  achieved  on  this  waveform.   So,  instead  of  trying 
to  detect  the  presence  of  the  signal  in  the  noise  using 
crosscorrelation  it  was  decided  to  make  simply  the  auto- 
correlation of  the  noiselike  waveform  at  the  output  of  the 
voltage  amplifier.   The  use  of  two  voltage  amplifiers  in 
series  became  necessary  in  order  to  overcome  a  well  defined 
input  threshold  level  since  only  after  this  signal  level  is 
achieved  the  equipment  can  perform  properly. 

A  non-expected  waveform  was  displayed  in  the  oscillo- 
scope after  the  autocorrelation  operation.   Rather  than 
being  associated  with  the  acoustic  waveform  expected,  it 
was  simply  the  result  of  the  autocorrelation  of  a  60  Hz 
power  line  interference  corrupted  with  random  noise.   It 
could  be  observed  that  the  same  autocorrelation  waveform 
was  exactly  obtained  whether  the  radar  was  transmitting  or 
not. 

As  a  conclusion  of  these  experiments  and  according 
to  the  theoretical  calculations,  such  a  low  input  electro- 
magnetic energy  is  not  enough  to  give  detectable  results 
after  suffering  this  extremely  weak  conversion  process. 
2.   Transient  Surface  Heating  Experiment 

In  the  second  group  of  experiments  performed  the 
fresh  water  was  assumed  an  elastic  medium.   If  results  were 
eventually  obtained  by  using  a  high  power  pulsed  laser 
they  could  be  attributed  mainly  to  a  thermally  shocked 
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interface  rather  than  to  the  simultaneous  but  much  weaker 
radiation  pressure  effects. 

Two  types  of  lasers  were  used  in  this  part  of  the 
experimental  work  as  sources  of  electromagnetic  energy. 

The  first  one  was  the  Nitrogen  pulsed  laser  existing 
in  Bullard  Hall  Room  217.   This  device  is  capable  of 
delivering  in  normal  operating  conditions  as  much  as  80 
kw  peak  power.   Its  pulse  repetition  rate  can  be  varied 
but  20  pulses  per  second  was  the  one  used.   The  pulse 
width  is  about  10  ysec  and  the  wavelength  is  0.3371  y. 

The  detection  apparatus  was  again  the  hydrophone, 
voltage  amplifier  and  oscilloscope. 

Unfortunately  the  output  light  of  this  source  has 
not  the  shape  of  a  narrow  beam;  instead  it  diverges  too 
much  and  in  order  to  focus  as  much  light  as  possible  several 
lenses  had  to  be  used.   A  mirror  was  also  used  to  deflect 
to  a  vertical  orientation  the  horizontal  output  beam.   As 
a  result  of  all  these  effects  and  characteristics  the  input 
energy  at  the  surface  of  the  water  was  not  enough  to  generate 
a  detectable  signal  at  the  oscilloscope  display.   Due  to 
improper  shielding  of  the  laser,  electromagnetic  interfer- 
ence at  the  pulse  repetition  frequency  was  picked  up  also 
random  noise  and  power  line  interference  were  observed. 

No  conversion  effect  was  then  detected  through 
this  scheme. 

The  second  laser  used  was  a  single  pulse  laser 
existing  in  the  Laser  Laboratory,  Building  223.   This  is  a 
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high  energy  Neodymium  Glass  Rod  laser  operating  at  the 
wavelength  1.06  y  and  its  output  light  is  highly  direc- 
tional.  The  beam  was  deflected  to  the  vertical  orientation 
using  a  low  absorption  prism.   It  is  possible  to  vary  its 
output  energy  and  values  up  to  12  Joules  can  be  safely 
obtained.   However  with  the  control  knobs  set  at  a  fixed 
position  it  is  observed  that  the  output  energy  do  not  stay 
constant  from  shot  to  shot.   The  extent  of  this  variation 
can  be  appreciated  next  in  Fig.  6. 

Through  the  use  of  a  beam  splitter  a  portion  of  the 
output  beam  was  diverged  toward  a  photo  detector  whose 
output  voltage  signal  was  used  to  measure  the  energy  of 
each  output  pulse.   Another  photo  detector  excited  by  the 
visible  light  of  the  laser  flash  tube  was  used  for  scope 
triggering  purposes.   This  triggering  action  took  place 
about  0.1  msec  earlier  than  the  beginning  of  the  laser 
pulse.   So,  there  is  a  time  delay  of  about  0.1  msec  between 
the  starting  of  the  sweep  and  the  beginning  of  the  laser 
pulse  that  should  be  considered  when  interpreting  Fig.  6. 
Two  TEKTRONIX  oscilloscopes  were  used:   a  memo-scope  type 
564B  to  measure  the  energy  of  each  pulse  and  a  545A  with 
single  sweep  capability  to  observe  and  photograph  the  de- 
tected signal  with  a  polaroid  TEKTRONIX  oscilloscope  camera 
system.   A  sketch  of  the  experimental  setup  is  shown  in 
Fig.  5. 
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FIG.  5  —  Experimental  Setup 

The  experiment  was  performed  with  sea  water  to 

decrease  the  penetration  depth  of  the  laser  beam.   In  fact 

14 
at  2.83x10   Hz,  which  is  the  laser  frequency,  the  pene- 
tration depth  for  sea  water  with  productivity  a  =  4  mhos/m 
is  about  1.19  cm  while  for  fresh  water  it  is  orders  of 
magnitude  greater.   Thus,  the  use  of  sea  water  would  guaran- 
tee a  total  absorption  of  the  electromagnetic  energy  within 
the  limited  system  used. 

The  hydrophone  was  not  exactly  in  the  same  vertical 
of  the  laser  beam  to  prevent  the  possibility  of  damage 
caused  by  the  intense  radiation. 

To  make  sure  that  the  signal  obtained  after  each 
shot  was  not  the  result  of  electromagnetic  interference, 
a  flat  piece  of  wood  was  placed  over  the  container  to  prevent 
the  radiation  from  reaching  the  surface  of  the  water. 
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(a)   6.6J,  no  input 


(b)   6.9J,  free,  5  cm 


(c)   6.6J,  free,  5  cm 


(d)   6.8J,  constrained,  5  cm 


(e)   6.3J,  free,  10  cm 


(f)   7.1J,  constrained,  10  cm 


(g)   5.7J,  free,  17  cm 


(h)   6.1J,  free,  17  cm 


FIG.  6  —  The  output  energy,  water  surface  condition  and 
hydrophone  depth  are  shown  for  each  shot. 
Horizontal  scale  is  0.1  msec/div  for  all  pictures 
Vertical  scale  is  5v/div  except  for  picture  (h) 
where  it  is  2v/div. 


Picture  6 (a)  was  taken  under  these  conditions  and  as  it 
can  be  seen  no  signal  was  recorded;  this  shows  clearly 
that  the  signals  obtained  before  came  really  from  a  surface 
heating  conversion  process. 

Fig.  6  presents  the  results  of  a  series  of  shots 
for .various  surface  conditions  and  hydrophone  depths .   The 
voltage  amplifier  was  kept  at  20  db  for  all  of  them.   As 
mentioned  before  although  the  laser  controls  were  at  a 
fixed  position  the  output  energy  was  not  constant  and 
values  ranging  from  5.7  to  7.1  Joules  were  obtained. 

From  pictures  6(b)  to  (f)  it  can  be  seen  that  the 
time  delay  between  the  starting  of  the  laser  pulse  and  the 
beginning  of  the  detected  signal  is  about  the  same  whether 
the  hydrophone  was  at  5  or  10  cm  depth  and  much  bigger 
than  the  time  that  sound  in  water  takes  to  travel  5  or  10 
cm.   Obviously  this  is  due  to  the  facts  that  heating  is 
not  instantaneous  and  that  the  energy  is  not  only  abosrbed 
at  the  surface  that  causes  the  thermal  effect  and  the  sound 
production  to  take  place  but  over  a  considerable  region,  the 
extent  of  which  is  determined  by  the  penetration  depth. 

Pictures  (g)  and  (h)  were  taken  with  the  hydrophone 
much  farther  from  the  beam  vertical  than  in  all  the  others 
and  at  17  cm  from  the  small  surface  area  illuminated.   Here 
the  effect  of  a  greater  distance  from  the  shocked  region 
is  apparent  since  the  time  delay  is  now  greater. 

Pictures  (d)  and  (f)  were  taken  under  constrained 
surface  conditions.   As  predicted  by  the  theory  the  signal 
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obtained  had  a  higher  amplitude;  it  also  happened  however 
that  the  energy  in  both  cases  was  higher  than  in  the  free 
surface  case  and  so  it  is  not  clear  in  which  extent  the 
increase  in  amplitude  was  determined  by  the  surface  constraint 

It  should  be  noted  that  a  circular  glass  plate  of  about 

2 

8  0  cm  was  used  touching  the  water  to  constrain  its  surface 

which  means  that  just  a  small  percentage  of  it  was  constrained 

Clear  is  the  fact  that  in  general  the  shape  of  the 
waveforms  obtained  is  not  close  to  the  theoretical  predictions 
made  before.   It  is  not  difficult  to  accept  such  a  different 
behavior  since  the  assumed  theoretical  conditions  were  not 
fulfilled  in  this  experiment,  namely  the  observation  point 
was  not  in  the  far  field  and  the  heating  was  not  at  or  near 
the  surface.   Moreover,  it  is  known  that  the  envelope  of 
the  output  pulse  of  this  laser  system  is  not  a  neat  rectangle 
but  exhibits  instead  two  peaks  and  has  a  shape  and  duration 
varying  from  pulse  to  pulse.   Also  the  relatively  long  plastic 
rod  supporting  the  hydrophone  may  easily  get  into  vibration 
when  subjected  to  the  acoustic  pressure  and  acoustical 
reflections  from  the  walls  are  likely  to  occur. 

To  appreciate  the  pressure  field  peak  amplitude 
consider  for  example  fig.  6(c)  where  a  peak  voltage  amplitude 
of  about  10  V  was  obtained.   Taking  into  account  the  20  db 
of  signal  amplification  and  using  eq.  (91)  it  can  be 
computed  that  the  pressure  amplitude  of  the  produced 
acoustical  field  is  3.16  bars  or  3.12  atmospheres. 
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B.   THEORETICAL  COMPUTATIONS 

Theoretical  calculations  related  to  the  radiation 
process  are  now  made  to  give  a  further  insight  on  the 
problem  and  to  allow  a  better  understanding  of  the  diffi- 
culty in  getting  experimental  results  at  low  input  power 
levels. 

Assume  that  the  medium  under  consideration  is  sea  water 
with  conductivity  a  =  4  mhos/m,  relative  permittivity  e  =  80 
and  relative  permeability  y   =  1. 

The  reflection  coefficient  R  for  a  plane  electromagnetic 
wave  incident  normally  on  the  surface  of  the  water  is  given 
by 


R  -  Er  _  [J  "  % 


where  E   and  E .  are  respectively  the  reflected  and  incident 
E  field  strengths,  n   =  /y  /e    =  3170,      is  the  intrinsic 
impedance  of  free  space  and 


n  =       -a«H 


CT  +  JU)£  £ 

J   o  r 


is  the  intrinsic  impedance  of  the  water.   When  the  frequency 

of  the  radiation  is   f  =  9.3  GHZ   then 

q 
oo  =  27Tf  =  58.43x10   rad/sec, 

n  =  42.02   |  2?76   , 

R  =  0.799   |  -0?623   . 
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This  means  that  only  20.1%  of  the  E  field  will  penetrate 
the  water. 

If  the  incident  electromagnetic  energy  is  from  a  radar 
with  output  peak  power  P   =  100  W  and  if  a  horn  antenna 
with  gain  G.  =  15.8  db  =  38.02  is  at  r  =  1  m  above  the 
surface  of  the  water,  then  the  power  density  at  the  interface 
is 

PtGt  2        Ei 

P.    =  -=-£  =  302.55   W/rn      =   — 
1       a     2  n 

4-rrr  'o 

so  E?  =  n   x  302.55   or  E.  =  337.73  V/m.   Then 
1     o  1 

337.73  x  0.799  =  67.88  V/m 

will  penetrate  the  water. 

The  attenuation  constant  of  the  water  is  given  by 


1/2 


Mei     A    x    ,_o,2 


a    =    oj[  J^M  x/1  +    (t^)      -   1)1  =      84.142   Np/m 

^    V         COG 


The  penetration  depth  <5  =  1/a  is  then  equal  to  1.18  8  cm. 
So,  as  far  as  absorption  of  energy  is  concerned,  the  medium 
may  practically  be  considered  as  extending  from  z  =  0  to 
z  =  0°.   Thus  the  components  of  the  radiation  pressure  force 
are  in  this  case  the  force  given  by  eq.  (4  6)  due  to  energy 
dissipation  of  20.1%  of  the  incident  energy  and  the  force 
given  by  eq.  (48)  due  to  the  reflection  at  the  interface  of 
79.9%  of  that  energy. 
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From  eq.  (46)  the  average  force  due  to  dissipation  of 
the  transmitted  energy  is 

E2 

rT        V        XO  a 

F,  =  — s-  x  cos  6 

vo    2M 


where  |n|  =  42.02  ft,  e  =  2? 76  and  E   =67.8  8  V/m;  v  is  given 

by  eq.  (24)  and  v  =  l//yF~  . 

2 

Computing  the  values  of  v  and  v   it  turns  out  that 


v  =  33.48  x  10   m/sec, 

v2  =  11.23  x  1014  . 
o 


The  average  force  is  then 


Fx  =  163.195  x  10"8  N/m2 

=  163.195  x  10"7  dyne/cm2 
=  163.195  x  10~7  ybars. 


With  this  value  for  the  average  pressure  field  and  knowing 
that  the  hydrophone  sensitivity  level  is  ML  =  -110  db 
relative  to  lV/ybar,  it  is  possible  using  eq.  (91)  to  compute 
what  would  be  the  hydrophone  output  voltage 

-110  =  20  log  V/f^ 
log  V/F1   =  -5.5 

V/F1   =  3.16  x  10"6 

V  =  0.516  x  10~10  Volts. 
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From  eq.  (48)  the  average  force  due  to  energy  reflection 
at  the  interface  is  given  by 


™  2s    i 


where  the  factor   /i  u     is  now  1  since  incidence  and 

r  r 

reflection  of  energy  are  taking  place  in  free  space, 

o 

c  =  3«  x  10   m/sec  and 

E2  2 

-  xo        269.85  QC    c__   TT/    2 

S   "  —=    2x377      =    96*577   W/m 


So 


F,   =  —  =    64.385    x   10"8   N/m2 
3  c 


=  64.385  x  10    ybars 


Following  the  same  procedure  as  before  the  output  voltage 
would  be 


V  =  F3  x  10~6  =  0.203  x  10  10  Volts. 


The  values  just  obtained  for  the  output  voltages  due 
to  both  force  components  are  so  small  that  they  cannot  be 
easily  seen  in  ordinary  oscilloscopes.   More  sophisticated 
detection  schemes  could  lead  to  positive  results  but  as 
mentioned  before  even  through  a  correlation  process,  it  was 
not  possible  to  detect  the  presence  of  the  signal. 
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Electromagnetic  energy  has  to  be  increased  orders  of 
magnitude  if  radiation  pressure  forces  are  to  be  easily 
displayed  by  detection  equipment. 
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IV.   CONCLUSIONS 

According  to  theoretical  predictions  and  experimental 
results,  in  order  to  get  significant  acoustic  power, 
extremely  high  electromagnetic  energy  has  to  be  used.   The 
reason  for  this  is  the  low  efficiency  that  characterizes 
this  type  of  energy  conversion. 

In  practice,  when  the  conversion  process  takes  place, 
both  radiation  pressure  and  transient  surface  heating  forces 
will  be  developed  in  the  medium. 

Thermal  forces  are  by  far  more  important  since  they 
are  orders  of  magnitude  larger  than  radiation  pressure 
forces.   The  present  status  of  laser  technology  has  spurred 
interest  in  the  possible  applications  of  this  conversion 
process.   In  the  future,  it  is  likely  that  there  will  be 
ocean  applications  of  highly  directive  thermoacoustic 
arrays  excited  by  modulated  lasers. 

Radiation  pressure  forces,  although  much  smaller  in 
magnitude  and  requiring  careful  design  of  detection  schemes, 
higher  power  sources  and  highly  sensitive  devices,  have  been 
proposed  to  be  used  in  several  applications. 

A.  Ashkin  [Ref .  17]  observed  that  micron-sized  particles 
in  liquids  and  gas  may  be  accelerated  and  trapped  in  a  laser 
beam.   The  experiments  were  performed  taking  special  care 
to  avoid  as  much  as  possible  thermal  effects  from  the  laser 
light  so  that  the  observed  effects  were  due  primarily  to 
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radiation  pressure  forces.   Chemical  and  laser  pumping 
applications  were  suggseted  on  the  basis  of  this  effect. 

H.E.M.  Barlow  points  out  in  his  paper  [Ref .  4]  the 
possibilities  of  the  application  of  radiation  pressure 
forces  in  power  measuring  instruments  and  in  orientation 
behavior  in  satellites. 

Much  work  has  to  be  done  yet  in  the  field  of  the 
conversion  of  electromagnetic  to  acoustic  energy  until  all 
useful  applications  are  thoroughly  explored. 

As  a  final  comment  it  should  be  noted  that  the  Laser 
Laboratory  at  NPGS  seems  to  have  the  appropriate  material 
available  to  encourage  further  experimental  work  on  the 
production  of  acoustic  waves  in  media  through  transient 
surface  heating. 
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APPENDIX  A 
DETAILED  THEORETICAL  CALCULATIONS 

1.   Penetration  depth  (6)  for  sea  water 

a   =   4  mhos/m 

e   =  80 
r 

-12 

£  =e  £   =  708.32  x  10    F/m 
o  r 

ur  =  1 

a.   For  the  radar  frequency   f  =  9.3  GHZ 

9         „„    „„  c 


CO 


=   2-jrf  =    2      x   9.3   x  10      =   58.43   x   10      rad/sec 


ue    =   u    e    e      =   4tt    x   10"7   x   8.854   x   10~12    x   80 
K  o  o   r 

=  8.901  x  10"16 


a    _  A _  4 


w£    58.43  x  109  x  708.32  x  10~12    41.3871 


=  0.0966484 


"-'JflJ1*  <^)2   "   l) 


=  58.43  x  109   n/4.4505  x  10~16  (  /l  +  0.009341    -  1) 


58.43  x  109   /"4.4505  x  10"16  x  0.0046596 


=  58.43  x  109   /  0.020738  x  10~16 
=  84.142  Np/m 

<S  =  i  =  0.01188  m  =  1.188  cm 

a 

83 


b.   For  the  laser  frequency 


g 

f  =  T  =   3  X  10 c    =  2.83  x  1014  HZ 

A    1.06  x  10"b 


a)  =   2Trf  =   17.78   x   10        rad/sec 


a 


"£    17.78  x  1014  x  708.32  x  lo"12    1259392.96 


=    3.1761   x   10"6 


14    a;     \\nr.        ,«-l6  .„-12 


a   =    17.78    x   10    ^    /   4.4505    x    10  x    5    x    10 

=    83.873   Np/m 


5   =  -  =    0.01192   m  =   1.192   cm 

a 


2.   Reflection  coefficient  (R)  for  sea  water  for  the  radar 
frequency 


e    n  -  n 
R  =  _£  =  J 2. 

Ei   n  +  nQ 


9 
oj  =  58.43  x  10   rad/sec 
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juy     m     I  J158.43  x  109  x  4tt  x  10~7 

4  +  J58.43  x  109  x  708.32  x  10"12 


J90 


o 


J73425.3035  /  73425.3035    e 

4    +    ^41'3871        ~J    41.5799e^84?48 


=  42.0224  ej2*76  =  41.974  +  J2.0235 


r    41.974  +  J2.0235  -  377  -335.026  +  J2.0235 

41.974  +  J2.0235  +  377  '  418.974  +  J2.0235 

-._  _  9  j0?3461 

=  0.799  r-0?623 


418.979e^°2767 


3.   Radiation  pressure  force  F,  due  to  energy  dissipation 

E2Q 

Fl.=  ^T  x  ifnT  C°S  9 
'  vo 


00 
V  =  y- 


ye  [ 


A   +  (a/\ie)2      +      1] 


8.901  x  10"16  [J   1  +  (0.0966484)2   +   l] 


8.901  x  10  16  x  2.00466 


=  0.3348  x  10   m/sec 


v  2  =  —  = rr  =  0.11235  x  1016 

°     ^    8.901  x  10"16 

-  =  0.3348  x  108    x   67. 882    x  Q^9QQ   =    ^^    x   1Q-8  N/m2 
0.11235  x  10      2  x  42.02 
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